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AEROTHERMODYNAMIC ENVIRONMENT OP A TITAN 
AEROCAPTURE VEHICLE 


By 

S. N. Tiwari 1 and H. Chow 2 
SUMMARY 

The extent of convective and radiative heating for a Titan entry vehi- 
cle is investigated. The flow in the shock layer ia assisted to be axisym- 
metric, steady, viscous, and compressible. It is further assumed that the 
gas is in chemical and local thermodynamic equilibrium and tangent slab 
approximation is used fcr the radiative transport. The effect of slip 
boundary conditions on the body surface and at the shock wave are included 
in the analysis of high-altitude entry conditions. 

The implicit finit^-dif ference technique is used to solve the 
viscous shock- layer equations for a 45-degree sphere cone at zero angle 
of attack. Different compositions for the Titan's N 2 + CH 4 atmosphere are 
assumed, and results are obtained for the entry conditions specified by the 
Jet Propulsion Laboratory. The results indicate that the heating rate, in 
general, increases with increasing N 2 concentration. Both convective and 
radiative heating increase with increasing initial entry velocity. The 
radiative heating increases, but the convective heating decreases with 
increasing body nose radius. The amount of CN concentration in the shock- 
layer gas determines the extent of radiative heating to the body. Radiative 
heating will be important for free-stream gas composition with N 2 concen- 
tration between 50Z and 90Z. For the atmospheric compositions of 99. 5X N 2 ♦ 
0.5Z CH 4 and 98Z N 2 ♦ 2Z CH 4 , the radiative heating hear the stagnation 
region is insignificant in comparison to the convective heating. The 
results indicate that the effect of the slip conditions is important when 
the altitudes are higher than 402.595 km. Therefore, both the body and 
shock slip conditions should be included in analyzing the aerothermal 
environment of the Titan aerocapture vehicle at higher entry altitudes. 

^Bainent Professor, Department of Mechanical Engineering and Mechanics, 

Old Dominion University, Norfolk, Virginia 23508. 

2 Graduate Research Assistant, Department of Mechanical Engineering and 
Mechanics, Old Dominion University, Norfolk, Virginia 23508. 
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1 . INTRODUCTION 


The importance of aerob raking and aerocapture for planetary missions 
has been emphasized in the recent years. It has been pointed out that aero* 
braking for circularizing orbits and aerocapture could more than double 
science payload on some planetary spacecraft and make possible new missions, 
such as a Saturn orb iter dual probe mission, where the probes would enter 
the atmospheres of both Saturn and its satellite Titan* 

The aerobraking technique uses the aerodynamic drag of the spacecraft 
during successive passes through the upper atmosphere to circularize a high* 
ly elliptical orbit. The aerocapture concept, on the other hand, uses the 
aerodynamic drag to place the spacecraft in a closed planetary orbit from a 
hyperbolic flyby trajectory in a single atmospheric entry pass. It is 
accomplished through an aerodynamical ly controlled atmospheric entry during 
which the vehicle's in-plane lift-t<^drag ratio is varied to maintain a 
constant drag. The aerocapture net only offers significant gains in payload 
and choice of orbits, but also significantly decreases interplanetary cruise 
time; and this concept completely eliminates the fuel-costly retropropulsion 
module for planetary orbiter mission. 

An aerocapture mission is possible for any atmosphere-bearing celestial 
body. The feasibility of using aerocapture vehicles has been emphasized 
recently for both inner and outer planetary missions (refs. 1-6). Origin- 
ally, the aerocapture study was undertaken for a Mars sample return mission 
(refs. 3,6). The aerocapture missions under present consideration are the 
Mars surface sample return (MSSR), Saturn orbiter dual probe (S02P), and 
Titan orbiter (TO) missions. 

For missions to outer planets, use of the aerocapture concept in a 
convenient atmosphere-bearing satellite of the target planet has been 
emphasized. It has been proposed to use the atmosphere of Titan for braking 
intc a Saturn orbit (ref. 4). The use of Titan's atmosphere would minimize 
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Che entry speed requirement for aerocapture and this, in turn, would mini- 
mize the thermal protection requirements of the aerocapture vehicle. The 
Titan's aerocapture concept (for Saturn orbital mission) is expected to cut 
the interplanetary cruise travel time to Saturn from 8 to 3.5 years. A 
Titan orbiter mission using anything other than aerocapture is presently 
impractical (ref. 4). For Titan's aerocapture mission, the need for high- 
performance entry vehicle geometries and high-performance thermal protec- 
tion systems has been stressed (refs. 4,5). In partial support of this 
need, it is essential to provide a complete analysis of the aerothermo- 
dynamic environment of the Titan aerocapture vehicle. 

The optimum lift/drag (L/D) ratio required for the aerocapture control 
accuracy is 1.0 to 2.0. The combination of high volumetric efficiency, low 
ballistic coefficient, and aerocapture control accuracy has led to choosing 
biconics as Che entry vehicle geometry for the aerocapture missions (ref. 

4). 

In order to investigate the aerothermodynamic environment of a Titan 
aerocapture vehicle, it is essential to know the composition of Titan's 
atmosphere. Prior to the Voyager 1 mission (November 1980), there was a 
controversy regarding Titan's atmospheric composition. The problem is still 
not completely resolved, but it is now evident (ref. 7) that Titan, the 
largest moon in the solar system, is wrapped essentially in a dense atmos- 
phere of nitrogen vapors (rather than methane, the best guess before Voyager 
1). Thus, a realistic composition for Titan's atmosphere would include a 
fairly high concentration of nitrogen. 

The main objective of this study is to determine the extent of convec- 
tive and radiative heating to the aerocapture vehicle under different entry 
conditions. This essentially can be accomplished by assessing the heating 
rate in stagnation and windward regions of an equivalent body. The equiva- 
lent body configuration considered for this study is a 45-degree sphere cone 
at zero angle of attack. Different compositions for the Titan's N 2 + CH^ 
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atmosphere have been asumed, and the study has been conducted for various 
entry trajectories suggested by the Jet Propulsion Laboratory (JPL)« Spe- 
cific obectives of this study, therefore, are as follows: 

1« For a given fre^stream atmospheric composition, determine the 
important chemical species in the shock- layer gas for different 
pressure and temperature conditions* 

2. Investigate the effect of the free-stream gas composition on the 
stagnation- point shock temperature and convective and radiative 
heating rates* 

3. Investigate the effect of different entry velocities on the stagna- 
tion-point shock temperature and convective and radiative heating 
rates* 

4* Determine the effect of body nose radius on the stagnation- point 
convective and radiative heating rates* 

5. Determine the variation of the shock temperature and enthalpy and 
convective and radiative heating rates along the body for different 
free-stream atmospheric compositions* 

6* Investigate the influence of CN concentration in the shock layer on 
the convective and radiative heating rates along the body* 

7. Investigate the effect of shock as well as body slip conditions on 
the entire shock- layer flow phenomena and determine the extent of 
convective and radiative heating rates under these conditions* 

Basic formulation of the entire problem is presented in Chapter 2, and 
boundary conditions are given in Chapter 3. The information on the thermo- 
dynamic and transport properties are given in chapter 4, and Chapter 5 dis- 
cusses the chemical compos itons* The radiative transport model for this 
study is described in Chapter 6. The physical conditions and data sources 
are given in Chapter 7* The method of solution is duscussed in Chapter 8, 
and all results are presented in Chapter 9* 
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2. BASIC FORMULATION 


As discussed in the Introduction, the aerocapture technique transfers 
the spacecraft into a closed, stable orbit from a hyperbolic flyby trajec- 
tory in a single pass (fig* 1)* This requires a high level of technology, 
but offers a significant gain in the payload and choice of orbits. The 
Titan aerocapture concept for the Saturn orbital mission is shown in figure 
2. For such missions, use of the biconics as the entry vehicle configura- 
tion (fig. 3) has been suggested (ref. 4). The preliminary assessment of 
the aerothermodynamic environment of an aerocapture vehicle can be made by 
investigating the flow field around an equivalent body. The equivalent body 
configuration considered for this study is a 45-degree sphere cone at zero 
angle of attack. 

The physical model and coordinate system considered for the equivalent 
body are shown in figure 4. The flow conditions for a radiating and reac- 
ting multicomponent gas mixture in the shock layer are considered axisymmet- 
ric, steady, viscous, and compressible. It is further assumed that the gas 
is in chemical and local thermodynamic equilibrium and the tangent slab 
approximation is used for the radiative transport. 

The conservation equations for a reacting multicomponent gas mixture 
can be found in the literature (refs. 8,9). The viscous shock-layer equa- 
tions that are valid uniformly throughout the shock-layer region are formu- 
lated in exactly the same manner as the viscous shock- layer equations for a 
one-component gas presented by Davis (ref. 10). In order to obtain the 
viscous shock- layer equations, the conservation equations are written in a 
boundary- layer coordinate system as shown in figure 4 and are nondimension- 
alized by variables which are of order one in the boundary layer. The same 
set of equations is then written in variables which are of order one in the 
inviscid region outside the boundary layer. Terms are kept in each set of 
equations up to second order in the inverse square root of Reynolds number. 
The two sets of equations are combined so that terms up to second order in 
both the inner and outer regions are retained. In this way, a set of 
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Figure 3. Titan aerocapture for Saturn orbit. 




equations uniformly valid to second order in the entire shock layer is 
obtained. The nondimens ional form of the viscous shock- layer equations that 
are applicable in the present case can be written as (refs. 11, 12): 


Continuity: 
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The terras used to nondiraensionalize the above equations are defined as 
u* » uV* v* » vV* 
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TV* 2 2 

T* • P* " P P * V* 




p* a p p* 


U* ■ UU* 


ref 


K* ■ Ku* .C* 
ref p 


c* • c c* 

p p p»* 


h* • hV* 


V* 

vrt • w ♦ * — 
1 1 R* 




J.M* , 
1 ref 




a* * aR^J 


n* » n& * 
N 


r* » rR* 
N 


K* - ICR* 


p* V* - 


Pr 


C*p* 
_P 

K* 


Le . . * p * C * 
P 
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JJ 
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( 2 . 8 ) 


The set of governing equation* presented above [eqs. (2.8)] has a 
hyperbolic-parabolic nature, where the hyperbolic nature come* from the 
normal equation. If the shock layer is assumed to be thin, then the normal 
momentum equation can be expressed as 


PW 2 < m 3p 

r 5n 


(2.9) 


When equation (2.3) is replaced with equation (2.9), then the resulting set 
of equations is parabolic. These equations can, therefore, be solved by 
using numerical methods similar to those used in solving boundary-layer 
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problems (ref*. 10,13). After an initial iteration using equation (2.9), 
the final flow-field solution is obtained by replacing equation (2.9) with 
equation (2.3); thus, the thin shock-layer approximation is removed. 


Since there are no nuclear reactions, the elemental mass fractions 
remain fixed and unchanged during chemical reactions. The relation between 
the elemental and species mass fractions is given by 




( 2 . 10 ) 


The elemental continuity equations for the elements can be obtained by 


multiplying equation 


* « 41 0 * 


ind sunning over i. The resulting 


elemental continuity equation is 




( 2 . 11 ) 


where 
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( 2 . 12 ) 


Use of the elemental mas* fraction reduces the number of equations to 
be solved. Therefore, equation (2.3) is replaced with (2.11) for equili- 
brium flow. 

The mas flux due to concentration gradients can be written as 




(2.13) 
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where 


Le. i - K 
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( 2 . 14 ) 
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The relative mass flux for the elements can be written as 



( 2 . 16 ) 
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and L is an arbitrary constant. For binary diffusion, equations (2.13) 
and (2.16) reduce, respectively, to 



(2.19) 


( 2 . 20 ) 


The heat transferred to the wall due to conduction and diffusion ia 
referred to here as the convective heat flux and is given by the relation 
(ref. 13): 



( 2 . 21 ) 


In this study, the Lewis and Prandtl numbers are taken to be 1.1 and 0.72, 
respectively. 

The convective heat transfer is also described by a dimensionless 
parameter called Stanton number. The Stanton number ia given by 


St » 



( 2 . 22 ) 


The skin friction coefficient for such flows is given by 



In order to solve the preceding set of governing equations, it is 
essential to specify appropriate boundary conditions at the body surface and 
at the shock. These are discussed in detail in the next chapter. 
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3. BOUNDARY CONDITIONS 


3.1. Introduction 

Specific boundary conditions used at the body surface and the bow shock 
are presented here. Since both the slip and no-slip conditions have been 
used in this study, they will be discussed separately in this chapter. 

3.2. No* SI ip Boundary Conditions 

At the body surface (wall), no velocity slip and no temperature jump 
are assumed. Consequently the velocities at the surface are 



V * 0 

(3.1) 


u - 0 

(3.2) 

The 

wall temperature for this study is specified as 



T * constant 
w 

(3.3) 

The 

surface total enthalpy is given as 



N 

H - y h.C. 
i-1 1 1 ' 

(3.4) 


The Rankine-Hugoniot relations are used to determine the flow proper- 
ties immediately behind the shock. The nondimens ional shock relations are 
as follows (refs. 10-13): 

Cont inuity: 


- sin a 


(3.5) 


Momentum: 


a" * cos a 
8 


P » + sin^a 1 

8 yj£ 



Energy: 


h 


s 





(3.6) 

(3.7) 


(3.8) 


State: 


p T R* 
s s 

P 8 = 

8 M* C* 


(3.9) 


Enthalpy: 

N 

h » T h.C. (3.10) 

8 i-i 1 1 

where a is shown in figure 4 and u" and v" are velocity components 

S 3 

expressed in a shock -oriented coordinate system. The transformations used 

to express u“ and v 11 in terms of the body-oriented coordinate system 
s s 

u 3 and v 8 are 

u 9 u" sin(a +0) + v M cos(a +0) (3*11) 

8 s s 


and 


v * -u M cos(a + 0) + v M sin(a + 0) 

8 3 s 


where the angle 0 is indicated in figure 4. 


( 3 . 12 ) 
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3.3* Slip Boundary Conditions 


In low Reynolds number hypersonic flows, such as high-altitude or low- 
density flows | the velocity and temperature of the wall are no longer the 
same as that of the gas immediately adjacent to the wall; these phenomena 
are referred to as the velocity slip and temperature jump, respectively. 

The slip flow boundary conditions have been derived by various investiga- 
tions (refs. 14-17). Shidlovskiy (ref. 14) has shown that at the body 
surface the velocity slip and temperature jump conditions are of the same 
order as the Knudsen number. The Knudsen number Kq is defined as the 
ratio of the particle's mean free path i and the characteristic dimension 
L of the body. These boundary conditions assume an impermeable surface and 
zero macroscopic velocity normal to the surface. They also assume that the 
mean free path i , although small, is large enough so that there is no 
interaction between incident and reflected molecules at the surface. Thus, 
for the transitional range, in order to be consistent with the Navier-Stokes 
equations of motion, a linear relation between the conditions at the wall 
and flow should be assumed. That this can be done is a semi-macroscopic 
argument which leads to simple expressions for the velocity slip and 
temperature jump as (refs* 14, 17): 




1/2 



(3.13) 

(3.14) 


T 


h 



where and A 2 are constants and are given by 


(3.15a) 


(3.15b) 
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,., 1/2 

A x - [(2 - O^/Ox] (-) 

2 

A 2 - 1(2 - o 2 )/o 2 ] (— )(-) 1/2 

8 2 

The terns cfj end o 2 are slip and thetmal accommodation coefficients, 
respectively, and are dependent on the nature of the surface and fluid. 
However, in actual flight conditions, both Oj and o 2 are expected to 
be 1. 

The boundary conditions used at the shock are the modified Rankine- 
Hugoniot or "shock slip" conditions, and these are written as (refs. 10, 
13): 



(3.16) 

(3.17) 

(3.18) 


+ sin 2 a - v" 2 ] (3.19) 

s 

u" sin(a + 8) + v" cos (o + 8) * u (3.20) 

8 8 3 

-u" cos(a + 8) + v" sin(a + 8) * v (3.21) 

8 8 8 

As mentioned above, slip boundary conditions are used in investigating 
the shock- layer flow phenomena at relatively high entry altitudes. 
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4. THERMODYNAMIC AND TRANSPORT PROPERTIES 


Thermodynamic properties for specific heat, enthalpy, and free energy 
and transport properties for viscosity, thermal conductivity, and diffusion 
coefficients are required for each species considered. Since the multicom- 
ponent gas mixtures are considered to be mixtures of thermally perfect gas- 
es, the thermodynamic and transport properties for each species are calcula- 
ted by using the local static temperature. The general expressions for 
total enthalpy, specific enthalpy, and specific heat at constant pressure 
are given, respectively, by 


H » h + 


(4.1) 


N 

h - l h.C. 
i«l 1 1 


(4.2) 


N 

C - l C.C . 
P iml i P.x 


(4.3) 


For each species , the values for the thermodynamic properties, as a 
function of temperature, are obtained by using polynomial curve fits. The 
following polynomial equations are used: 

Specific heat: 


C* 


R* 


~ “ a l + % T* + a 3 T* 2 + T*2 + ag T**» 


(4.4) 


Enthalpy: 


h* 

i 


a + *n + + ^ t * 3 + hl 

R*T* 1 2 3 4 5 T* 


(4.5) 
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\ 


L 


Free energy: 



R*T* 


ai(l “ logg T*) - 


2 3 

« 2 T* _ a 3 T» _ ai»T* 

2 6 12 
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T* 


a7 


(4.6) 


where F*0 is the free energy of species at one atmospheric pressure. 

The development of these curve fits and the values of polynomial constants 
a 1 to are given in table 1 and are available in reference Id. 


For the mixture, viscosity and thermal conductivity are obtained by 
using the semi-empirical formula of Wilke (ref. 8) as 



where 



(4.7) 


(4.8) 


(4.9) 


The general relations for the viscosity and the thermal conductivity are 
given as 

Hi * bi + b2 T* + b3 T* 2 (4.10) 
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(4.11) 


» Cj + C 2 T* 

The coefficients bp b 3 , Cp end C 2 for different species used in this 
study are given in table 2, where the value of T* is in degrees K. 
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5. CHEMICAL COMPOSITION 


Analyses of chemically reacting flows are usually simplified by assum- 
ing the chemical equilibrium behavior of the gas mixture. In this study, 
the chemical reactions are confined to a system of carbon, hydrogen, and 
nitrogen. The Aerotherm Chemical Equilibrium (ACE) computer program was 
used to determine various chemical species under different free-stream 
atmospheric compositions. 

At the initiation of this study, the atmospheric conditions of Titan 
were not defined clearly. Therefore, different atmospheric compostions were 
assumed for a parametric study. Voyager 1 data reveals that Titan ( s atmos- 
phere primarily consists of nitrogen molecules (ref. 7). Thus, a realistic 
case would be to assume a very high concentration of nitrogen in the free- 
stream gas mixture. However, to study the effect of free-stream gas compo- 
sition on heating of the entry vehicle, different gas compositions are as- 
sumed. 

The equilibria chemical composition is determined by using a free 
energy minimization analysis as developed in reference 19. As mentioned 
above, the ACE cooputer program was used to determine various chemical 
species for different pressure, temperature, and free-stream conditions. 

For initial study, 68 chemical species for the carbon-hydrogen-nitrogen 
system were included in the matrix of calculations for a given free-stream 
atmospheric compostion. The matrix was 

Pressure: 0.1, 0.5, 1.0 and 0.5 atm) 

Temperature: 2,000 K to 10,000 K in 500 K increments 

Composition: 90% N 2 + 10". Ch 4 , 502 N 2 + 50% C^ , 10% N 2 + 90% CH^ 

For different free-stream gas compositions, the variation in mole frac- 
tion of different species, as a function of temperature is illustrated in 
figures 5(a) to (d) for different pressures. There are about 20 chemical 
species shown in these figures. However, concentrations of some species are 
less than 0.05 percent for the range of temperature considered. Therefore, 
for this study. 17 chemical species (N 2 , N, N + , C 3 , C 2 1 C, C + , C 4 H, C 3 H, 
c 2 * h .> C 2 li * CN, ^2 * HCN, and E") were considered for the shock-layer 

gas mixture. 


0 2 4 6 8 10 X 10 

Temperature, K 

(a) 


Figure 5. Variation in mole fraction of different species for p * 
atm and 90% N? + 10% CH 4 . 













MOLf 

FRACTION 
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(c) 


Figure 5. (Continued.) 


30 




4 6 8 10 X 10 

Temperature, K 

(d) 


Figure 5. (Concluded.) 




6. RADIATION TRANSPORT MODEL 


An appropriate expression for the radiative flux, q r , is needed for 
the solution of the energy equation presented in Chapter 2. this requires a 
suitable transport nodal and a meaningful spectral model for variation of 
the absorption coefficient of the gas. 

In the present analysis the "tangent slab" assuaption for radiative 
transfer has been used. This implies that the radiative energ/ transfer 
along the body is negligible in comparison to that transferred in the direc- 
tion normal to the body. It should be noted the tangent slab approximation 
is used only for radiative transport and not for other flow variables. For 
a nonscattering medium and diffuse nonreflecting bounding surfaces, a one- 
dimensional expression for the spectral radiative flux is given by (refs. 

20 , 21 ): 


qrvUv) ■ 2ir {e v (B v (0)E3(x v ) - B v (t ov )E3 ( t ov - T v ) 1 



0 


B v U)E2(t v - t)dt 


T 

/° V B v (t)E2(t - t v )dt> 

T 


( 6 . 1 ) 


where 


y 

\ • / « v (y')dy' 

o 
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The quantities B y (0) and B y (T 0V ) represent the radiosit lee of the body 
surface and shock respectively. The expression of total radiative flux is 
given by 

<l r " r « rV < T v> dw (6 * 2) 

o 

In the shock layer, the rediative energy from the bow shock usually is 

neglectd in comparison to the energy absorbed and emitted by the gas layer. 

The expression for net rediative flux in the shock layer, therefore, is 
given by combining equations ( 6 . 1 ) and ( 6 . 2 ) as 

q r - 2 r h v (0)E 3 (T v ) ♦ / V B v (t)E 2 (T v - t)dt 

O 0 

- /° V B v (t)E 2 (t - T v )dt]dv (6.3) 

T 

r 


where q v (0) ■ e v vB v (T s ). 

In this equation, the first two terms on the right represent the radia- 
tive energy transfer towards the bow shock while the third term represents 
the energy transfer towards the body. Upon denoting these contributions by 
q r * and q r ~, equation (6.3) can be written as 



The radiative flux, q r is calculated with the radiative transport 
code RAD (ref. 22) which accounts for detailed nongray radiation absorption 
and emission processes. The chemical species considered for determining the 
radiative transport are H, M 2 , N*, IT, N^, H, H 2 , r, e", C, (T, C * , C 2 , 

Cj , and CN, 
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7. PHYSICAL CONDITIONS AND DATA SOURCE 


As mentioned earlier! the entry body considered for this study is a 45- 
degree sphere cone at a zero degree angle of attack (fig. 4). The body 
temperature is assumed to be 2,000 K and, for most cases, the body nose 
radius is taken to be 0.2 m. The free-stream atmospheric compositions are 

assumed a* 99.5 X N 2 ♦ 0.5X CH* , 98% N 2 + 2% CH 4 , 90X N 2 ♦ 10% CH 4 » 75% N 2 + 
25% 50% N 2 + 50% , 25% N 2 + 75% Cl^ , and 10% N 2 + 90% CH^. The high 

nitrogen concentration case., will be the realistic con*>ositions for the 
Titan's atmosphere. 

For th< Titan aerocapture mission, entry trajectories have been gener- 
ated by JPL. The altitude history for an aerocaputre vehicle is illustrated 
in figure 6 for two different (shallow and steep) entry angles. The entry 
trajectories and free-stream conditions used in this study are given in 
tables 3 to 9. 



Figure 6. The altitude history for an aerocapture vehicle. 



8* METHOD OF SOLUTION 


A numerical procedure for solving the viscous shock-layer equations for 
stagnation and downstream regions is given by Davis (ref. 10). Moss (ref. 
13) and Tiwari and Szema (ref. 12) applied this method of solution to react- 
ing multicomponent mixtures. A modified form of this procedure is used in 
this study to obtain solutions of the viscous shock-layer equations. In 
this method, a transformation is applied to the viscous shock-layer equa- 
tions in order to simplify the numerical computations. In this transforma- 
tion most of the variables are normalized with their local shock values; the 
transformed variables are (refs. 12, 13): 





P 

P 



u 

s 




K 

s 




T 

s 



V 


V 


V 

s 



8 


( 8 . 1 ) 
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The crane forma tions relacing Che differencial quantities are 

3 ( ) _ 3j_) _ n dn e 3 ( ) 

3x 3? n d£ 3n 

s 

3 ( ) _ 1 3 ( ) a 2 _ 1 3 2 ( ) 

3n n 3 3n n 3h 

8 S n 8 8 


( 8 . 2 ) 


(8.3) 


After the governing equations are written in the transformed variables, the 
resulting second-order partial differential equations can be expressed in 
the following form: 


3 2 W 
3n 2 


«i 


3W 

— + ♦ as 

3n 



0 


(8.4) 


The quantity W represents u in Che s-momentum equation, 1 in 
the enthalpy energy equation, in the elemental continuity equa- 

tion. The coefficients ai to at* to be used in this study are exactly 
the sane as given in references 12 and 13. 


8 -momentum, W * u: 


n k 

3 


n cos 6 
s 


1 3tt 

ai * + 

_3n 1 + nntc r + nn cos 6 

tr s s 


n p u n' — n p v 

s s s s pun _ 8 s s pv 


e 2 p (1 + n nic) p e 2 p p 
s s s 


(8.5a) 


«2 


<n , K 2 n 2 cos 9 n 2 < 

8 1 dp S s 


(1 + n ten) p 3n (l + n n<) 2 (r + n n cos 0 ) (1 + n n*c) 
s s 8 s 


P n 2 u‘ 
3 s 

e 2 p (1 + n nie) 
s a 


up 

u 


n 2 p v < 

8 3 3 £V 

e 2 p (1 + n n<) jT 
s s 


(8.5b) 
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(8.5c) 


P n 2 
~8 8 


3P . P 8^ n s n 3P 


e 2 p u (1 + n nic ) 


r 35 P 


8 8 


n 3n 

8 


a* » - 


pun 2 

8 8 8 PU 

e 2 p (1 + n nic) V 
8 s 


(8.5d) 


Energy (enthalpy, W * H: 


1 3p 1 3Pr 

a l*=— = +n s‘ 

U 3n Pr 3n 1 + n n< r + n n cos 6 


k cos 6 

+ 


p Pr Pr n n u nup 
8 s s ss -- 

■ - V pv 

e 2 p V 1 + n ntc 8 


(8.6a) 


H' 

8 

«2 * ^ — 
U 

8 


(8.6b) 


*3 


Pr Pr n 2 

3 8 1 3 \|) 


< cos 8 


* 


Prp v v 
*s s 3P 


p pH n3n 1 + n n< r + nn cos 0 e^p pH 3n 


s s 


8 8 


PrP 


r.s 


1 3 «r 


e 2 u H u n 3n 

S3 S 


+ <l D 


cos 9 


(8,6c) 


1 + n x\k) r n cos 0 
s s 


Pr n^p u Prup 

S 8 3 6 

e 2 u (1 + n n<) v 

a s 


(8.6d) 
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where 


n Pr 
8 8 


iT p . ^i u 2 uu , _ 3u 

_ 2. h . + — — - (Pr Fr - 1) — 

Tr i-1 1 3n Fr 3n 


V u^uu 2 
8 S 

1 + n n< 
8 


(8.6e) 


The preceding energy equation is for the thin shock-layer approximation. 
When equation (2.3) is used for the n-momentum equation, the following term 
must be added to equation (8.6c): 


Pr v n z p - — 
s s s s Prpv 


e 2 p H IT 
s s 


uu / , - n'nv.-\ 

8 /— , 3v _ s s 3v I 

VV s + v s — / 

1 + n nic \ 3£ n 3n / 


2 ” 

V*V , - 

s 3v 

n 3n 
s 


(8 . 6 f ) 


Elemental continuity, W * Cj^ : 


1 3 PL 

*1 + n 

~ 3n 
PL 1 


cos 9 


1 + nntc r + nn cos 9 
a s 


p v n v n p u n*pun 
s sss ^ ssss 


e 2 PL e 2 PL ( 1 + n nic) 

s 


(8.7a) 


«2 “ 0 
1 

“ — 

PL 


3 ~ ~ / k cos 9 \ 

— PM + n PM [ + j 

3n 8 \1 + n nic r + nn cos 9/ 

\ S 3 ' 


(8.7b) 


(8.7c) 
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i 


(8.7d) 


«4 


n 2 p u 
8 8 8 


c z (i + n n<) 

s 


PU 

PL 


where 


and 


U P. 

PL- JJl 


Pr 7r 
8 


★ 

NS M„ NS 3C. 

PM - I 6 .1 I Ab — 
i-1 * k-1 1K 3n 

M i *1 


for multicomponent diffusion and for binary diffusion: 


P PC 

PL - 

Pr Pr 
s 


PM - 0 

The remaining equations are written as follows: 
Continuity: 

— [ n (r + n n cos 0)p upu] + — [(r + n 
3£ 8 8 88 3n 

{(1 + n nic) p v pv - n'np u pu}1 « 0 
s s s s s s J 


s 


n cos 6)] 


(a.7e) 


(a.7f) 


(8.7g) 


(8.7h) 


(a. a) 


n-momentum: 


pu 

1 + n we 

s 




v — . — 
a pv 3v 

u n 3n 
8 8 


_ U 8 K — 2 . ^8 3P 

pu* + - 0 

v (1 ♦ n t\k) p u n v 3n 

8 8 8 8 8 8 


(8.9a) 

which becomes 


• 

_ n p u 2 < 

3P 8 8 8 — 2 

— ■ pu* 

3n p (l + n ntc) 

8 8 


(8.9b) 

if the thin shock~layer approximation is made 

. 


State: 



* 

F-pT L 
★ 

M s 


(8.10) 

The boundary conditions at the body surface (the 
conditions) in terms of transformed variables are as 

surface boundary 
follows: 

No slip: u * 0 


(8.11a) 

v « 0 


(8.11b) 

T • const. 


(8.11c) 

N N 

H - I h.C . / T h.C . 

i-1 1 1 w i-1 118 


(8. lid) 
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With slip: 




(8.12a) 

(8.12b) 

(8.12c) 

(8.12d) 


The conditions at the shock (i.e., the transformed shock conditions at n 
* 1) for slip or no slip cases are 


u“T«H*v-P«p*l 


(8.13) 


When downstream numerical solutions are required, it is necessary to 

have an accurate solution for the flow along the stagnation streamline. A 

truncated series which has the same form as that used by Kao (ref. 23) is 
1 

used to develop the stagnation streamline equations. The flow variables are 
expanded about the axis of symmetry with respect to nondimens ional distance 
£ near the stagnation streamline as 


p(5,n) * Pj(h) + P2(h)5 2 + ... 

(8.14a) 

u(£,n) * uj(nK ♦ ... 

(8.14b) 

v(£,n) * vj(n) + ... 

(8.14c) 

p(Cih) * p^n) + ... 

(8.14d) 

T(£,n) ■ Tjtn) + ... 

(8. 14e) 
(cont'd) 
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h(5,n) » hi (n) + 


u(5,n) « ui^) + . . 


k(£ ,n ) - Ki ) + ... 


c (5,n) - c (n) + ., 

p p.i 


c t (C,n) - c t x (n ) + ., 


(8.l4f) 

(8.14g) 

(8.14h) 

(8.14i) 

(8.14j) 
(concl 'd) 


The shock-standoff distance is expressed by 

n - n + n 5 2 + . . . 

* 1 .® 2 ® 


(8.15) 


Since 5 is small and the curvature k is approximately of order one 
in the stagnation region, it is logical to say that (see fig. 4): 


6-5 


(8 16) 


Now, since 6 * (ir/2) - 6, one may express 

/ 2rt -s 


o - 1 + 5 j - - 1| 

2 V * “i* 


(8.17) 


By using equations (8.15) to (8.17), the shock relations, equations 
(3.3) to (3.8), can be expressed in terms of expanded variables as 



(8.18) 


(8.19) 
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p. " Pi. ♦ P2.* 2 * ••• 


1 



( 8 . 20 ) 


( 8 . 21 ) 


Since equations (8.19) and (8.20) involve n 2 § , these terms cannot be 
determined from the stagnation solutions. Thus, a value of n 2 s ■ 0 is 
assumed to start the solution. This assumption is removed by iterating on 
the solution by using the previous shock-standof f distances to define n^. 

Along the stagnation streamline, the second-order differential equation 
is mitten as 


d 2 W dw 

+aj. — + a2 + 83 »0 

dn 2 dn 


( 8 . 22 ) 


The coefficients in equation (8.22) are defined as 
s -mo men turn, W ■ u: 
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(8.23a) 


(8.23b) 
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, - -"u^u h . 

'■‘“i. 1 * "l. n) “u "■ P 1. n i. *> 


(6.23c) 


Energy (enthalpy), W ■ 9; 


du dPr 

4l - 1 1 - 1 1 + 


2n n p Pr v p.Pr v 

le _ la la la la 1 11 


F^ dn Vr^ dn 1 ♦ n^ n 


e 2 u 


(6 1 24a) 


la 


a2 ■ 0 


(8.24b) 
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"1/1. Pr i i » 


U. H. F, n dn 1 + n n 
la la 1 la Is 
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(8.24c) 


Elemental continuity, W « : 


~ In npvpv 

1 dPL ^ ? la - ls ls Isj; 1 

~ dn 1 1 + nn, 

PL L !• e z PL 


(8.25a) 


a2 ■ 0 


1 

a 3 " ~ 
PL 


*U 2 ° 1<m 1 
dn 1+nn la J 


(8.25 b > 


(8.25c) 


The remaining equations are written as follows: 


Continuity: 


d_ 

dn 
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(i + 


n n ) 2 
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(8.26) 


45 


n -momentum: 


dF p dv 

1 l$ p ls 1 


(8*27) 


When the thin shock-lsyer approximation is made* the n-momentim equation 
becomes 


(8*28) 


The governing second-order partial differential equations are solved by 
employing an implicit finite-difference method* A variable grid spacing 
(fig. 7) ia uaed in Che n-direction so chat the grid spacing can be made 
small in Che region of large gradients. In Che figure, m is a sCaCion 
measured along Che body surface and n denotes Che station normal Co Che 
body surface. The derivatives are converted to finite-difference form by 
using Taylor's series expansions. Thus, unequal space central difference 
equations in the n-direction at point m, n can be written as 


3n 0 * + A V ”’ n+1 An n-l (An n-l+ An n> “ , “" 1 


An n ~ An „-1 

♦ n 01 w_ 


A V"n-l •*" 

32 wx 2 2 

3n 2 n An n (An n ♦ An n-1 ) W *’ n * 1 An^n^ “ ,n 


+ An n-l (An n + ^o-l > Wm *° +1 


(8.29a) 


(8.29b) 




(8.29c) 


W - u 

3W\ u s»,n nrl,n 

H* AC 


A typical finite-difference expansion of the standard differential equation 
is obtained by substituting the above equations in equation (8.4) as 

AW +BW +CW + D * 0 (8.30) 

n m, n+1 n m,n n m, n-1 n 


where 


B 

n 


2 + SjAri 

An + An 
n 


n-1 


n-1 



a, (An - An ,) 
1 \ n n-1/ 

An An . 
n n-1 



2 - a An 
1 n 

A Vi( An » * A Vi) 


(8.31a) 


(8.31b) 


(8.31c) 
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- a 




4 m-l,n 


nrl 


(8.3 Id) 


If it is assumed that 

W - E W + F (8.32) 

m,n n m, n+1 n 


or 


W 

m, 


n-1 


E W ♦ F 
n-1 m,n n-1 


(8.33) 
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Chen substituting (8.33) into equation (8.30) yields 


W -j ~ An w 

m,Q I B + C E °* n+1 

In n n-1 


-D -C F 
n n n-1 

B + C E , 
n n n-1 


(8*34) 


By comparing equations (8.32) and (8.33), one finds 


E 

n 


-An 


B - C E , 
n n n-1 


(8.35) 


F 

n 


-D - C F , 
n n n-1 


B 

n 


C E . 
n n-1 


(8.36) 


Now, since E. and F, are known from the boundary conditions, E and 
11 n 

F can be calculated from equations (8.35) and (8.36). The quantities 
n 

W n Q at point m,n can now be calculated from equation (8.32). 

The overall solution procedure starts with evaluation of the flow pro- 
perties immediately behind the shock by using the Rankine-Hugoniot rela- 
tions. With known shock and body surface conditions, the solutions are 
obtained first for the stagnation streamline. With this solution providing 
the initial conditions, the solution is marched downstream to the desired 
body location. Each of the second-order partial differential equations is 
integrated numerically by using the tridiagonal formalism of equation (8.4) 
and following the procedure described by equations (8.30) to (8.36). The 
first solution pass provides only an approximate flow-field solution. This 
is because, in the first solution, pass, the thin shock-layer form of the 
normal momentum equation is used, the stagnation streamline solution is 
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assumed to be independent of downstream influence, the term dn^/dC is equal 
to zero at each body station, and the shock angle a is assumed to be the 
same as the body angle 6. These assumptions are removed by making oue or 
more additional solution passes. 


The shock solution procedure at any location is identical for the first 
and subsequent solution passes. However the shock angle a is defined 
differently for the first and subsequent solution passes. For the first 
solution pass, a » 0. For subsequent solution, the shock angle is defined 
as 

a « 0 + tan 


1 + Kn 


(8.37) 


In the first solution pass, the viscous shock-layer equations are 
solved at any location m after obtaining the shock conditions. The 
converged solutions at station m - 1 are used as the initial guess tor the 
solutions at station m. The solution is then iterated -locally until 
convergence is achieved. For the stagnation streamline, guess values for 
dependent variables are used to start the solution. 


In the first local iteration, both dn^/d^ and 3w/€ are assumed to be 
zero. The energy equation is integrated numerically to obtain a new temper- 
ature. By using this temperature, new values of thermodynamic and transport 
properties are calculated. Next, the x-momentum equation is inegrating to 
find the u component of velocity. The continuity equation is used to 
obtain both the shock-standoff distance and the v component of veloc- 
ity. The pressure p is determined by integrating the normal momentum 
equation. The equation of state is used to determine the density. For 
example, the integration of the stagnation streamline continuity equation 
from 0 to n results in 


[(1 + n l8 n) 2 P ls v l8Pl ] v x - (-2n ls P l8 u l8 )A 


(6.38) 
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where 


n 

A • f (1 + n, n)p.u dn 
Is 1 1 

0 

This equation give the v-velocity component along the stagnation streamline. 
However, integration of the continuity equation from n ■ 0 to n ■ 1 results 
in 


Is 


Is Is 




(8.39) 


where 


1 

a * / p^udn, c 


l 

n, / p.u ndn 
Is 11 
o 


The shock-standoff distance can be obtained from the solution of equation 
(8.39) as 


Is 


- (2v , +2Bu, ) + [ ( 2v +2Bu ) 2 -4(v, +2Cu )v ] 1/2 
Is Is Is Is Is Is Is 

2( '’U +2C “1 S > 


(8.40) 


Integration of the downstream continuity equation from n * 0 to h results in 


3 r n — 

— LJ n (r+n n cos 6) p u pudn] 
o 4 am sm * * 


s s 


+ (r + n n cos e)I(l ♦ nn tc ) (p v pv) “ n 1 np u puj * 0 (8.41) 

sm sm s s sm s s 
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This can be expressed in terns of the difference equation as 


[(GG) - (GG) ,] 

m nr- X 

A? 


+ (FF) v + (EE) » 0 
m in 


(8.42) 


where 

(EE) * (r + n n cos 0)(1 + n n<)p v p 
m sm 8m s s 

(FF) * -(r + n n cos 8)n' np u pu 
m sm sm s s 

(GG) * f n n (r v n n cos 0)p a pudn 
m o sm sm 8 s 

Now, the v-velocity component at each point on the station m can be 
obtained from equation (8.42). 

For the downstream shock-standoff distance, integration of the 
continuity equation from n * 0 to n * 1 gives 

- — [n^ cos 0 p u j^pundn + n rp u f^pudn 
s s s ' s s s 4 

3£ o o 

* (r + n cos 0 ) [n P u - (1 + n <)9 v ] (8.43) 

a s s s a a a 

By defining, for station m 

D. * cos 0 p u f\>undn. D ■ rp u f^pudn 
1 s a 2 s a * 

o o 


and denoting the same relations by Dj an d 04 for station m- 1 , equation 
(8.43) can be expressed in terns of a difference equation as 


“V . 2 * “iV. - <D 3°. * \ 


r P u n* + cos Opun'n - r p v 
ss sm sssmsm as 


- r p v <n - cos 0 p v n - cos 0 p v <n 

ss sm ss sm 8 8 s m 


This can be expressed in a quadratic form as 


(II)n + (JJ)n + (KK) 
sm sm 


where 


I 1 * D. + cos 0 k p,v,4( 


JJ * D„ + r p v - cos fl p u o' 
2 s s s s s 


KK - - [D,(n ) + D, (n ) . + r p u n'4? - r P v 4§] 

3 a m-1 4 s mrl 8 s s 8 a 


Then, the shock-standoff distance at station m is obtained from equation 
(8.45) as 


n - (-(FF) ♦ ((JJ ) 2 - 4(II)(KK) ] 1/2 } [2(H )]” 1 
am 


The flow diagrams for computation procedure are shown in figures 8 (a) 
to 8 (d). 
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Figure 8(c). Flow chart for subroutine energy solution 
procedure. 






Figure 8(d) . 


Flow chart for subroutine momentum 
solution procedure. 









9. RESULTS AND DISCUSSION 


The entry body considered for thie study is s 45-degree sphere cone st 
zero degree angle of attack. The body surface temperature is taken to be 
uniform at 2,000 K and the body nose radius is 0.2 m. The entry trajector- 
ies and free-stream conditions are given in tables 3 to 9. Results have 
been obtained to investigate the effects of different gas composition, entry 
velocity and body nose radius on the stagnation- point convective and radia- 
tive heating. Specific results were obtained to determine the extent of 
convective and radiative heating along the body for free-stream gas composi- 
tion of 90% N 2 ♦ 10% CH 4 , 98% « 2 ♦ 2% CH„, and 99.5% N 2 ♦ 0.5% CH 4 . Jhe 
results for Trajectories I to VI with different free-stream gas compositions 
are given in tables 10 to 18. For the slip boundary conditions, some impor- 
tant results are presented in this section to show the effects of both the 
body and shock slips on the convective nesting; these results are given in 
tables 19 to 24. 

For Trajectory I, the effects of free-stream gas composition on the 
shock temperature are illustrated in figures 9 and 10. The results show 
that the shock temperature as well as the temperature in the shock layer 
increases with increasing N 2 concentration. This is because N 2 provides 
less energy accomodation in comparison to CH The stagnation shock temper- 
atures are relatively higher for early entry time (fig. 10); this, however, 
would be expected because of relatively higher free-stream velocities. The 
results of figure 9 show that the temperature gradient in the shock layer is 
restricted essentially in the regions near the body surface for all free- 
stream gas compositions. 

The effects of gas composition on the stagnation-point convective and 
radiative heating for Trajectory I are illustrated in figures 11 and 12. 

The convective heating is seen to increase with increasing N 2 concentration 
(fig. 11), and peak heating occurs at an entry time of about 70 s. This is 
a direct consequence of the variation in the shock temperature. The situa- 
tion, however, is not the same with respect to the radiative heating (figs. 


58 


Figure 9. Effect of gas composition on temperature distribution along 
the stagnation streamline. Trajectory I (time ■ 78 s). 
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Figure 10. Effect of gas composition on stagnation-point shock 
temperature. Trajectory I. 
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Figure 12. (Concluded.) 
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12(a) and (b)], i.e.) the radiative heating does not necessarily increase 
with increasing N 2 concentration. This is because, for a given set of con- 
ditions, the radiative transfer strongly depends on the presence of absorb- 
ing-emitting species in the gas mixture. It is also evident from figures 
12(a) and (b) that the peak radiative heating occurs at different entry 
times for different free-stream gas compositions. For N 2 concentrations 
between 50Z and 90Z; the maximum heating is noted for 75Z N 2 concentration. 
The free-stream compositions of 99. 5Z N 2 + 0.5Z Ch 4 and 98Z N 2 + 2Z CH 4 are 
considered to be the realistic compositions for Titan's atmosphere. The 
results for stagnation-point shock temperature and convective and radiative 
heatings are illustrated in figure 13 for Trajectory I and for 99. 5Z N 2 + 
0.5Z CH 4 . It is noted that for this case, the radiative heating is negligi- 
ble as compared to the convective heating. The radiative heating is not 
more than seven percent of the total heating. 

The effect of entry velocity on the stagnation- point shock temperature 
and convective and radiative heating rates are illustrated in figures 14 to 
16. For the free-stream gas composition of 90Z N 2 + 10Z CH 4 , the results 
presented in figures 14(a), 15(a) and 16(a) show that the shock temperature 
and heating rates, in general, increase with increasing entry velocity for a 
fixed entry altitude (time). It is seen that the extent of convective heat- 
ing is considerably higher than the radiative heating for all cases. The 
results also show a similar trend for the gas composition of 99 .5Z N 2 + 0.5 
CH 4 [figs. 14(b), 15(b), and 16(b)]. For this gas composition, the radia- 
tive heating is negligible in comparison to the convective heating. One 
exception to this, however, is noted from the results presented in figure 
16(b). The radiative heating rate for an entry velocity of 13 km/s is con- 
siderably higher than for other velocities. Thus, for the entry speed of 13 
km/s (and for an entry time between 30 and 60 s), it is possible to have 
physical conditions in the shock layer to produce a higher concentration of 
radiating species. 

For the free-stream atmospheric composition of 90Z N 2 + 10Z CH 4 , the 
results for stagnation- point convective and radiative heating are shown in 
figure 17 (and also in table 14 for Trajectories I and II), The results 
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Figure 13 . 
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Figure 14(a). Effect of entry velocity on stagnation-point shock 
temperature, 90% N 2 + 10% CHi». 
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Figure 14(b). Effect of entry velocity on stagnation-point shock 
temperature , 99.5% N 2 +0.5% CHi,. 
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Figure 15(a). Effect of entry velocity on stagnation-point convectiv; 
heating, 90% N 2 + 10% CH 4 . 
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Figure 15(b). Effect of entry velocity on stagnation- point convective 
heating, 99. 5% N 2 + 0.5% CH*. 
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Figure 17. Verurioo of etagnetloo-polnt convective and radiative 
heating for Trajectories I end II. 901 Nj + 101 CH». 




•how that the extent of both convective and radiative heating ie consider- 
ably higher for Trajectory I (a steeper entry angle trajectory) than for 
Trajectory II. This, however, is expected since the rate of viscous dis- 
sipation will be higher for the steeper trajectory, resulting in a relative- 
ly higher shock temperature. 

For the atmospheric composition of 99. 5% N 2 ♦ 0.51 CH 4 , the variation 
of at agnation* point convective and radiative heating with body nose radius 
is given in figure 18 for Trajectory I and an entry time of 78 s (2 • 50.9 
km). Although the extent of radiative heating is small, it is seen to in- 
crease with increasing nose radius. The convective heating rate, however, 
is seen to decrease with increasing nose radius. For a given set of entry 
conditions, the shock- standoff distance generally increases with increasing 
nose radius (ref. 4). This, in turn, results in different temperature, 
pressure, and species distributions in the shock layer. A combination of 
these changes influences the trend exhibited. 

The results of heating rate along the body are illustrated in figures 
19 and 20 for Trajectory I and for the entry conditions at 78 s. Variations 
in shock temperature and heating rates are shown in figure 19 for the atmos- 
pheric composition of 99.52 N 2 * 0.5% CH 4 . The results show that both con- 
vective and radiative heating esaen ^lly follow the trend of the shock 
temperature from the stagnation point to the tangency point (at about S/R* * 
0.8). Beyond this point, the convective heating continues the same trend, 
but the radiative heating is seen to increase with the body location. This 
is because the pressure and temperature conditions near such location are 
conducive for production of the radiating CN species over a large portion of 
the shock-layer thickness (see fig. 5), and also because the optical 
thickness of the shock-layer gas is relatively higher in the downstream 
regions. The variation in heating rates along the body is illustrated in 
figures 20(a) and 20(b) for the cases with and without CN concentration in 
the shock-layer gas. The results show that while the presence of CN has 
little influence on the convective heating, the radiative heating is incre- 
ased considerably by its presence. It is important to note that after the 
tangency point, the rate of radiative heating in the presence of CN is sign- 
ificantly higher than the convective heating for the free-stream composition 
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Figure 18. Effect of body nose radius on stagnation-point convective 
and radiative heating for 99.5 X N 2 + 0.5Z CHu, t • 78 a. 
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Figure 20(a). Influence of CN on convective and radiative heating 
along the body for t « 78 s and 90% N 2 + 10% CH 4 . 



of 90% N 2 + 10% CH 4 [(fig* 20(a)]* The same trend is seen in figure 20(b) 
for the atmospheric composition of 99. 5Z N 2 + 0.5Z CH 4 , but the extent of 
radiative heating is considerably small* 

For the free-stream atmospheric composition of 99. 5Z N 2 + 0*5Z CH 4 , 
variations in important results with distance along the body surface are 
illustrated in figures 21 to 23 for Trajectory I and for critical entry 
times (altitudes). The results for shock density and shock- standoff dis- 
tance presented in figure 21 show that t for a given entry altitude, the 
shock- standoff distance increases as density decreases. The shock-standoff 
distance is seen to decrease with increasing altitude; this is because high- 
er free-stream velocities are associated with higher altitudes (see table 
1). The results for shock temperature and enthalpy presented in figure 22 
show that both decrease along the body until the tangency point, and they 
remain essentially constant beyond that point. Because of higher free- 
stremn velocities, the shock temperature and enthalpy are greater for higher 
altitudes. The variation in heating rates is shown in figure 23. As dis- 
cussed earlier, the peak convective heating occurs for entry conditions at t 
* 70 8 (2 * 70.4 km) and peak radiative heating at t * 78 s (z * 50.9 km). 

These results clearly show that the radiative heating is not important in 
the stagnation region if Titan's atmospheric composition is considered to be 
99. 5Z N 2 + 0.5% CH 4 . 

Trajectory VI is the latest trajectory specified for the Titan mission, 
and there appears to be a general agreement to consider the atmospheric 
composition as 98% N 2 + 2% CH 4 . Results for this case are shown in figures 
24 and 25. The results for stagnation- point shock temperature, enthalpy, 
and convective and radiative heating rates are shown in figure 24 for dif- 
ferent entry times. The results show that the extent of radiative heating 
for. this trajectory is small compared to the convective heating. It is 
noted that the radiative peak heating occurs at an entry time of 70 s (Z * 
204.570 km), whereas the convective peak heating (20.276 MW/m 2 ) occurs at 73 
3 (Z * 196.349 km). The variations in shock temperature, shock density, 
shock- standoff distance, and heating rates along the body are shown in 
figure 25 for an entry time of 73 s. These results exhibit essentially the 
same trend as noted in figures 19 and 21 for Trajectory I with 99.5% N 2 + 
0.5% CH 4 . 
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Figure 21. Variation of shock density and shock-standoff distance with 
body coordinate for 99*5% N 2 +0.5% CH 4 . 



Figure 22. Variation of shock temperature and enthalpy with body 
coordinate for 99.5% N 2 +0.5% CHi». 
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Figure 24. Variation of stagnation-point shock temperature, enthalpy 
and convective and radiative heating for Trajectory VI, 
98% N 2 + 2% CH 4 . 
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Thus, for atmospheric compositions with very high N 2 concentration, the 
radiative heating is not important in the stagnation region. 

The extent of convective and radiative heating over the entire length 
of the aerocapture vehicle is shown in figure 26 for the free- stream gas 
composition of 90% N 2 ♦ 10% CH 4 and 98% N 2 ♦ 2% CH 4 . For the free- stream 
composition of 90% N 2 ♦ 10% CH 4 , the results clearly show that, while the 
convective heating rate continues to decrease in the downstream region, the 
radiative heating rate is considerably higher in this region. As discussed 
before, the reason for this trend is the combined influence of shock- temper- 
ature density and pressure variations in this region and the relatively 
higher optical thicknesses of the radiating shock layer. A similar trend in 
heating rates is noted also for the gas composition of 98% N 2 + 2% CH 4 , but 
the extent of radiative heating is found to be relatively lower. However, 
it is important to note that the radiative heating approaches the convective 
heating in the downstream region. 

To investigate the effects of body and shock slip conditions on the 
entire shock- layer flow phenomena, the results were obtained for the recent- 
ly specified Trajectory IV (Table 6), with a free-stream atmospheric compos- 
ition of 98% N 2 + 2% CH 4 . Since chemical equilibrium is assumed and the 
thickness effect is of higher order, the concentration slip and thickness 
effects were neglected in this study. The results were obtained specifical- 
ly for the higher altitude entry conditions where the influence of slip 
conditions was anticipated. Some important results of this investigation 
are presented here. In discussion of these results (and in figures), the 
word ,, slip ,r implies both the body and shock slip conditions. Results are 
presented first for the velocity slip and temperature jump at the body 
surface. Following this, results are presented for the properties immedi- 
ately behind the shock. Next, the effects of slip conditions on the temper- 
ature distribution in the shock layer and on the convective heating along 
the body are discussed. The results are then presented for selected entry 
altitudes to show the separate effects of body and shock slip conditions on 
the convective heating along the body. Finally the results are presented 
for the temperature distribution and convective heating for very high entry 
altitudes . 


84 


. MW/m 






The variation in the surface-slip velocity is illustrated in figure 27 
as a function of the entry time. The results clearly illustrate that the 
condition of no slip is not satisfied at higher altitudes (earlier entry 
time). The magnitude of velocity slip (u wls ) is expressed as a percent 
of the velocity just behind the shock. It is evident from the figure that 
about 7 percent velocity slip occurs at the entry time of 20 s (2 ■ 465.115 
km), and only 0.12 percent at the entry time of 73 a (i ■ 196.349 km). 

The temperature jump and enthalpy change at the body surface are shown 
in figure 28 for different entry altitudes. The body surface temperature 
was taken to be 2,000 K. A temperature jump of about 1.2 times the surface 
temperature (i.e, AT * 2,500 K) is seen along the b surface near the 
stagnation point for entry altitude Z * 465.115 km. At lower altitudes, 
however, the temperature jump is seen to be relatively smaller. For exam- 
ple, at Z * 196.349 km the temperature jump is about 150 K. A similar trend 
is noted for enthalpy change at the body surface. 

Figures 29 to 32 show the velocity slip, temperature jump, enthalpy 
change, and density change just behind the shock. The results illustrated 
in figure 29 show that both the u and v velocity components are influ- 
enced by the slip conditions. It is evident from figures 30 to 32 that, 
when the altitude is lower than 402.595 km, the effects of slip conditions 
are not important. However, a significant temperature jump is noted at Z * 
465.115 km (see fig. 30). Since both the temperature and velocity compo- 
nents decrease just behind the shock, the slip conditions result in a 
decrease in enthalpy and an increase in density; these are clearly evident 
from the results in figures 31 and 32. 

The temperature distribution in the shock layer, along the stagnation 
streamline, is shown in figvre 33 for different altitude entry conditions. 

It is evident from the figure that, when the altitude is lower than 402.595 
km, the effect of slip conditions is not important. 

The effects of slip conditions on the convective heating along the body 
are shown in figure 34 for different entry conditions. The effects are seen 
to be lower for lower altitudes. It is important to note that at lower 
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Figure 27. 


Velocity slip at the body surface as a function of entry 
time (altitude) on the stagnation point. 
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Figure 30. Temperature variation just behind the shock wave as a 
function of £ coordinate for different entry altitudes 
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Figure 34. Variation of convective heating along the .body surface for 
different entry altitudes. 



altitudes the slip conditions result in an increase in the convective heat-* 
ing, whereas the reverse is true for the higher altitudes. For example, 
about a 1 percent increase in convective heating is noted for Z ■ 196.349 km 
and about 6 percent increase for Z * 241.638 km; however, a reduction of 
about 48 percent is observed for entry conditions at Z ■ 465.113 km. 

Separate effects of the body and shock slips on the 'onvective heating 
are shown in figures 35 and 36 along with the slip (body as well as shock 
slip) and no slip solutions. As would be expected, the results obtained by 
considering only the body or shock slip fall, in general, between the re- 
sults of slip and no-slip conditions. The effects of slip, of course, are 
higher for higher altitude entry conditions. The results clearly indicate 
that both the body and shock slips are equally important in influencing the 
extent of convective heating to the body. 

In order to assess the effets of slip conditions on the convective 
heating for very high entry altitudes, the results were obtained by consid- 
ering only the body-slip conditions because of the computational conveni- 
ence. The temperature distribution along the stagnation streamline is il~ 
lastrated in figure 37 for entry conditions at Z * 497.656 km and Z * 

531,004 km. The results show that, while the no-slip tem>erature distribu- 
tion is essentially the same for both altitudes, the body-slip temperature 
distributions are entirely different. The wall temperature jump and convec- 
tive heating variation along the body are illustrated in figure 38. These 
are seen to be influenced greatly by the slip-body conditions. For entry 
conditions at Z * 531.004 km, the results show a temperature jump of about 
150 percent and a decrease in convective heating by about 30 percent. 

The results for slip conditions clearly indicate that both the body and 
shock slip conditions should be included in analyzing the aerothermal envi- 
ronment of the Titan's aerocapture vehicle at higher entry altitudes. How- 
ever, during most of the heating pulse (where the heating is significant 
conpared with peak heating), this study indicates that accurate results can 
be obtained without including slip boundary conditions while using the as- 
sumption of equilibria flow. 
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Figure 36. Variation of convective heating along the body for different 
slip conditions at Z » 402.595 km and Z « 465.115 km. 



Figure 37. Temperature profile In the shock layer at stagnation point 
with body slip condition for very high altitudes. 




10. CONCLUSIONS 


The main objective of thie study was to assess the extent of convective 
and radiative heating that would be experienced by an aerocapture vehicle in 
a Titan mission. Different compositions for Titan's atmosphere were asstmd 
and results were obtained for the entry trajectories specified by JPL. The 
influences of slip boundary conditions (both at the shock and the body) were 
investigated for important cases. Specific results were obtained for free- 
stream atmoapheric compositions of 90X N 2 + 10X CH 4 , 99. 5Z N 2 ♦ 0.5X CH V , 
and 98Z N 2 ♦ 2X CH V . 

Results show that both the convective and radiative heating rates are 
quite sensitive to the gas composition used. The convective heating in- 
creases significantly as the N 2 concentration increases. However, this, in 
general, is not the case with regaid to the radiative heating. The radia- 
tive heating is negligible for the shallow entry (y * -25*) condition re- 
gardless of the free-stream gas composition. But, for the steepest entry 
angle (Y ■ -45*), the radiative heating will be important only if the free- 
stream gas is assumed to contain N 2 concentrations between 50Z and 90Z. For 
the gas composition of 90Z N, , the radiative heating is important in the 
stagnation region (as well as in the downstream region) with the peak radia- 
tive heating rate being 30 percent of the corresponding convective heating 
rate (about 13 MW/m 2 ). For the free-stream gas composition of 99. 5X N 2 ♦ 
0.5Z CH 4 , the radiative heating, in the stagnation region, is negligible 
(less than 6Z) in comparison to the convective heating for all cases consid- 
ered. For this gas composition, the peak convective heating is found to be 
about 15 MW/m 2 . The amount of CN concentration in the shock* layer gas 
determines the extent of the radiative heating. For a given free-stream gas 
composition, the radiative heating downstream of the stagnation region in- 
creases due to an increase in the CN concentration and the optical thickness 
of the shock layer. 

Other results obtained in this study show thst higher initial entry 
speeds produce higher shock temperature which, in turn, results in higher 
hesting rate. Results for the gas composition with 99. 5Z N 2 indicate that, 
while the convective heating decreases, the radiative heating increases with 
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increasing body nose rsdius. Specific results obtained for a more recent 
trajectory (Trajectory VI with free- stream gas composition of 98Z N 2 ♦ 2Z 
CK 4 ) indicate that radiative heating becomes comparable to the convective 
heating in the far downstream region of the aerocapture vehicle. 

Results of slip conditions (for the recent trajectory with free-stream 
gas composition of 98Z 1» 2 ♦ 2% CH 4 ) clearly indicate that both body and 
shock slip are important in influencing the shock-layer flow phenomena for 
high- altitude conditions. As such, these should be considered in determin- 
ing the extent of heating rates to the aerocapture vehicle at higher entry 
altitudes. 

For further study, it is suggested to consider the influence of chemi- 
cal as well as radiative nonequilibrium in analyzing the aerothermal envir- 
onment of the aerocapture vehicle. At this time, it might be advisable also 
to include the effects of thickness and concentration slip. However, during 
most of the heating pulse, this study indicates that accurate results can be 
obtained without including slip boundary conditions while using the assump- 
tion of equilibrium flow. i . 
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Table 2. Viscosity and thermal conductivity constants. 


SPECIES 

bl 

b 2 

b3 


ci 

C 2 

N 

.253000E-05 

. 220600E-07 

-.373700E-12 

* 

.128100E-04 

•859300E-08 


.9 OOOCE-05 

•161300E-07 

-. 19160CE-12 

* 

.6 54 00 (£-05 

.645700E-08 

N+ 

0 . 

. 50000 0E-08 

-.100000E-12 

* 

.260000E-03 

0 . 

c 

.19970CE-04 

♦177200E-07 

-.337800E-12 

it 

.25060CE-04 

.747900E-08 

<52 
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•139300E-07 

-.257500E-12 
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°3 
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Table 3. Altitude and free-stream conditions: Trajectory 1 



(L/D 

- 1-2, Y - 

-45*, 8 - 

800 kg/n? # 

U E - 10 km/s). 

TIME 

ALTITUDE 

P “3 



v. 

MACH 

(s) 

(km) 

(g/cm ) 

(mb) 

(K) 

(km/s) 

NO. 

40 

230.965 

0. 1265E-6 

0.10355 

159.10 

9.929 

29.13 

50 

169.824 

0.3473E-6 

0.24562 

139.44 

9.803 

30.72 

60 

114.238 

0. 1022E-5 

0.64627 

119.17 

9.431 

31.90 

70 

70.409 

0.3197E-5 

1.5836 

103.16 

8.448 

30.78 

78 

50.922 

0. 5157E-5 

2.5901 

96.28 

7.185 

27.10 

90 

48.539 

0.5512E-5 

2.7449 

95.56 

5.502 

20.83 

100 

60.654 

0.3944E-5 

1.9738 

99.26 

4.727 

17.56 

150 

95.613 

0.1604E-5 

0.89555 

112.46 

3.406 

11.89 

220 

173.831 

0. 3728E-6 

0.23159 

140.84 

3.015 

9.40 


Table 4. 


Altitude and free-stream conditions: Trajectory II 

(L/D - 1.2, y * -25*, 6 - 600 kg/m 2 , U - 10 km/s). 

£ 


TIME 

(s) 

ALTITUDE 

(km) 

P “3 

(g/cm ) 

(mb) 

T. 

00 

v. 

(km/s) 

MACH 

NO, 

50 

321.264 

0.3952E-7 

0.3249E-1 

160.00 

9.953 

29.11 

100 

213.217 

0.1654E-6 

0 . 1335E+0 

155.60 

9.689 

28.74 

110 

198.471 

0. 2063E-6 

0. 1638E+0 

150.30 

9.583 

28.91 

120 

185.306 

0.2694E-6 

0.1967E+0 

145.10 

9.451 

29.03 

130 

173.617 

0.3288E-6 

0. 2323E+0 

140. 76 

9.286 

28.96 

140 

163.335 

0.3812E-6 

0.2683E+0 

137.16 

9.100 

28.75 

150 

154.339 

0.468 7E-6 

0. 3139E+0 

134.02 

8.689 

28.41 

170 

141.502 

0.6712E-6 

0.3909E+0 

129.53 

8.358 

27.18 

175 

139.412 

0. 7052E-6 

0.4053E+0 

128.79 

8.210 

26.77 

180 

137.610 

0.7213E-6 

0.4215E+0 

128.16 

8.062 

26.36 

200 

130.832 

0.7855E-6 

0.4875E+0 

125.79 

7.495 

24.73 


no 


7 

j 


Table 5. Altitude and free-stream conditions: Trajectory III 



(L/D 

- 1.2, Y - 

-45*, 8 • 800 kg/m 2 , 

Ug ■ 6 km/a) 

l. 

TIME 

ALTITUDE 

P "3 

0m 

T. 

v. 

MACH 

(s) 

(lun) 

(g/cm ) 

(mb) 

(K) 

(km/a) 

HO. 

90 

144.674 

0.6142E-6 

0.3719 

130.00 

5.847 

18.93 

100 

111.891 

0.1076E-5 

0.6732 

118.70 

5.681 

19.29 

110 

83.284 

0.2349E-5 

1.1359 

108.10 

5.379 

19.14 

120 

61.655 

0.3839E-5 

1.9338 

99.60 

4.872 

18.06 

129 

49.822 

0.5318E-5 

2.6615 

95.90 

4.327 

16.35 

140 

45.041 

0.6079E-5 

2.9723 

94.50 

3.690 

14.04 

150 

47.699 

0.5643E-5 

2.7995 

95.30 

3.245 

12.30 

160 

54.496 

0.4666E-5 

2.3577 

97.30 

2.934 

11.01 

170 

63.709 

0.3693E-5 

1.8516 

*100.50 

2.724 

10.05 
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Table 6. Altitude and free-stream conditions: Trajectory IV 

(L/D - 1.2, Y - “45*. 0 - 800 kg/m?, U £ - 8 km/a). 


TIME 

(•) 

ALTITUDE 

(km) 

P “3 

(g/cm ) 

(mb) 

T. 

(K) 

V. 

(km/s) 

MACH 

NO. 

60 

181.218 

0.2927E-6 

0.2070 

143.48 

7.883 

24.35 

70 

134.830 

0.7469E-6 

0.4464 

127.19 

7.709 

25.30 

80 

93.931 

0. 1689E-5 

0.9283 

111.87 

7.329 

25.63 

90 

63.296 

0.3726E-5 

1.868 

100.32 

6.515 

24.07 

97 

50.882 

0.5162E-5 

2.592 

96.26 

5.783 

21.81 

110 

46.913 

0.5768E-5 

2.850 

95.07 

4.554 

17.28 

\20 

54.841 

0.4621E-5 

2.335 

97.17 

3.961 

14.84 

130 

67.354 

0.3413E-5 

1.706 

101.94 

3.603 

13.20 

150 

98.351 

0.1473E-5 

0.8421 

113.42 

3.240 

11.25 
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Table 7. Altitude and free-stream conditions: Trajectory V 

(L/D - 1.2, y - ”45* , 0 - 800 kg/in 2 , U £ - 13 ka/a). 


TIME 

(s) 

ALTITUDE 

(tan) 

P **3 

(g/cto ) 

P. 

(ab) 

T. 

<K) 

v. 

(ka/s) 

MACH 

NO. 

30 

237.512 

0.11438*6 

0.0937 

159.75 

12. 905 

37.77 

40 

158.551 

0.4165E-6 

0.2887 

135.49 

12.672 

40.28 

50 

90.769 

0. 1863E-5 

0.9899 

110.77 

11.777 

41.40 

60 

51.374 

0.5092E-5 

2.560 

96.41 

9.361 

35.28 

61 

49.676 

0.5339E-5 

2.671 

95.90 

9.084 

34.33 

70 

49.557 

0.5357E-5 

2.678 

95.86 

7.047 

26.64 

80 

67.533 

0.3401E-5 

1.698 

102,01 

5.912 

21.66 

100 

101.803 

0.1345E-5 

0.7892 

114.72 

5.083 

17.56 

120 

124.659 

0.8488E-6 

0.5381 

123.63 

4.762 

15.84 


Table 8. Altitude and free~streaa conditions: Trajectory VI 

(L/D - 1.4, y - -36* , 8 - 800 kg/m 2 , U £ - 12 km/a). 


TIME 

(a) 

ALTITUDE 

(km) 

P **3 

(g/cm ) 

0- 

(mb) 

x. 

00 

v. 

(km/e) 

MACH 

NO. 

10 

531.004 

0.1282E-08 

0.846E-03 

177 

11.994 

44.09 

15 

497.656 

0.2366E-08 

0.185E-02 

177 ' 

11.996 

44.10 

20 

465.115 

0.4441E-08 

0.432E-02 

177 

11.997 

44.11 

30 

402.595 

0.1489E-07 

0.904E-02 

177 

11.995 

44.10 

60 

241.838 

0.4396E-06 

0.4279E 

177 

11.664 

42.88 

70 

204.570 

0.1019E-05 

0.7567 

177 

11.096 

40.79 

73 

196.349 

0. 1227E-05 

0.8292 

177 

10.843 

39.86 

77 

187.956 

0.1483E-05 

0.9033 

177 

10.459 

38.45 

81 

182.652 

0. 1672E-05 

0.9501 

177 

10.043 

36.92 

90 

181.525 

0.171 5E-05 

0.9600 

177 

9.144 

33.62 

100 

192.414 

0.1341E-05 

0.8640 

177 

8.405 

30.90 

110 

200.836 

0.1109E-05 

0.7897 

177 

7.902 

29.05 

120 

204.913 

0. 1012E-05 

0.7537 

177 

7.510 

27.61 

130 

205.759 

0.9926E-06 

0.7462 

177 

7.175 

26.37 
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Table 9. Free- *t ream thermodynamic values for different 
gae compositions. 


MOLE ] 

FRACTION 

MIXTURE 

MOLECULAR 

HEIGHT 

(M) 

MIXTURE 

SPECIFIC 

HEAT 

(ft*/s - R) 

MIXTURE 

ENTHALPY 

(Btu/lbm) 

10% 

♦ 

90Z 

CH^ 

17.238 

11816.0 

1933.22 

25Z N* 

♦ 

75Z 

CHh 

19.035 

10884.6 

1622.04 

50Z 

♦ 

50Z 

CH4 

22.030 

9326.5 

1103.42 

75Z N* 

♦ 

25Z 

CH4 

25.025 

7771.2 

584.80 

90Z N* 


10Z 

CH 4 

26.822 

6836.5 

273.62 

98Z N/ 

♦ 

2Z CH4 

27.780 

6338.2 

107.66 

99. 5Z N* 

♦ 0 

>5Z CH4 

27.960 

6235.2 

76.54 
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Table 10. Stegnacion results (sphere cone, ■ 0.2 m, T^ ■ 2,000 K.) : 
atmosphere - 99.5X N 2 ♦ 0.5X Cft,, Trajectory 1. 


(a) 


TIME 

P- 

P . 

P,/P* 

T . 

“a 

q c,m 

q r,w 

(r) 

(kg/m 3 ) 

(atm) 


00 

(cm) 

(MW/m 2 ) 

(MW/m 2 ) 

40 

0. 12656-3 

0.1163 

18.13 

8508 

0.7963 

5.783 

0.0044 

50 

0.3473E-3 

0.3112 

17.93 

8498 

0.8241 

8.558 

0.0155 

60 

0.1022E-2 

0.8473 

17.76 

8262 

0.8381 

12.265 

0.0576 

70 

0.3197E-2 

2.1194 

16.75 

7902 

0.8927 

14.818 

0.2391 

78 

0.5157E-2 

2.4560 

15.08 

7354 

1.0011 

11.188 

0.3136 

90 

0.5512E-2 

1.5121 

12.10 

6412 

1.2341 

4.688 

0.2249 

100 

0.39448-2 

0.7888 

10.48 

5813 

1.4109 

2.290 

0.1387 


1 


Table 10. (Concluded. ) 


(b) 


TIME 

(e) 

(kJ/kg) 

h . 

(kJ/kg) 

SC 

Re 

40 

2. 043E+03 

48.799E+03 

0.9797E-0) 

0.106X+04 

50 

2.045E+03 

47.629E+03 

0.5497E-01 

0.2867E+04 

60 

2.04X+03 

44.082E+03 

0.301X-01 

0. 8317E+04 

70 

2.045E+03 

35.275E+03 

0.1640E-01 

0.2438E+05 

78 

2. 04.9! +03 

25.445E+03 

0. 128CE-01 

0.378X+05 

90 

2.045E+03 

14.799E+03 

0.1200E-01 

0. 3679B*05 

100 

2.045E+03 

10.873E+03 

0. 1374E-01 

0.251X+O5 


Table 11. Stagnation results (sphere cone, ■ 0.2 m, • 2,000 K) : 
atmosphere - 99. 52 N 2 + 0.52 CU 4 , Trajectory III. 


(a) 


TIME 

(s) 

P« 

(kg/m? ) 

P s 

(atm) 

P s /P» 

T s 

00 

n s 

(cm) 

^c,w 

^r, w 
(MW/m 2 ) 

90 

0.6135E-03 

0.1926 

13.96 

5983 

1.0951 

2.035 

0.0110 

100 

0.1075E-02 

0.3174 

13.31 

6056 

1.1407 

2.416 

0.0232 

110 

0. 2345E-02 

0.6164 

12.22 

6103 

1.2286 

2.895 

0.0669 

120 

0.3855E-02 

0.8218 

10.84 

5908 

1 .3689 

2.545 

0.1323 

129 

0.5312E-02 

0.8805 

9.45 

5546 

1.5510 

1.860 

0.2084 

140 

0.6072E-02 

0.7185 

8.10 

4823 

.1.7890 

1.009 

0.217 ‘ 

150 

0.5637E-02 

0.5113 

7.56 

4004 

1.9090 

0.528 

0.1709 

160 

0 .4660E-02 

0.3448 

7.37 

3358 

1 . 9600 

0.275 

0.0889 
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Table 11. (Concluded.) 


(b) 


TIME 

(s) 

h v 

(kJ/kg) 

h 8 

(kJ/kg) 

SC 

Re 

90 

2.045E+03 

16. 764E+03 

0.382CE-01 

0.453X4-04 

100 

2.045E+03 

15.840E+03 

0.2844E-01 

0.7693E+04 

110 

2. 045E+03 

14. 1 6 IE f 03 

0. 1876E-01 

0.159X4-05 

120 

2.045E+03 

11.550E+03 

0.1406E-01 

0 .2488E+05 

129 

2.045E+03 

9.045E+03 

0.1135E-01 

0.3262E4-05 

140 

2.045E+03 

6.512E+03 

0.9842E-02 

0.3618E4-05 

150 

2. 045E+03 

4. 982E+03 

0.9514E-02 

0. 345X4-05 

160 

2.045E+03 

4.038E+03 

0.9694E-02 

0 .2984E4-05 


Table 12. Stagnation results (sphere cone, ■ 0.2 m, ■ 2,000 K) : 
atmosphere - 99 . 5 % M 2 + 0.52 CH 4 , Trajectory IV. 


(a) 


TIME 

(s) 

(kg/m* ) 

P s 

(atm) 

P 8 /P. 

T . 

(K) 

“s 

(cm) 

^c.w 

(MW/m?) 

q r,* 
(MW/m 2 ) 

60 

0. 2923E-03 

0. 1696 

18.07 

6658 

0.8557 

4.039 

0.0047 

70 

0 . 7460E-02 

0.4219 

17.26 

6942 

0.8921 

6.032 

0.0191 

SO 

0. 1687E-02 

0.8396 

16.05 

7004 

0.9536 

7.233 

0.0589 

90 

0.3722E-02 

1 .4520 

14.27 

6866 

1.0610 

7.118 

0.1549 

97 

0.5156E-02 

1.5700 

12.71 

6561 

1. 1805 

5.527 

0.2061 

110 

0.57O0E-02 

1 .0640 

9.95 

5759 

1 .4780 

2.385 

0.2291 

120 

0.4616E-02 

0.6349 

8.67 

5147 

1.6800 

1.205 

0.1618 

130 

0.3410E-02 

0.3844 

8.03 

4640 

1 .8060 

0.667 

0.1037 


120 


Table 12. (Concluded.) 


(b) 


TIME 

(s) 

h „ 

(kJ/kg) 

h 8 

(kJ/kg) 

St 

Be 

60 

2.042+03 

30. 707E+03 

0.608CE-01 

0.2472+04 

70 

2.045E+03 

30.027E+03 

0.3680E-01 

0.6063E+04 

80 

2.0454+03 

26.496E+03 

0.2970E-01 

0.1302+05 

90 

2.045E+03 

20.873E+03 

0.1540E-01 

0.2681E+05 

97 

2.042+03 

16.374E+03 

0.1279E-01 

0.3516E+05 

110 

2.045E+03 

10.051E+03 

0.1117E-01 

0.3588E+05 

120 

2.042+03 

7.534E+03 

0. 1 176E-01 

0. 2784E+05 

130 

2.045E+03 

6.193E+03 

0.1275E-01 

0.2050E+05 


Table 13. Stagnation results (sphere cone, ■ 0.2 m, » 2,000 K) 
atmosphere - 99. 5% N 2 + 0.5% CH^, Trajectory V. 


(a) 


TIME 

(s) 

Poo 

(kg/m 3 ) 

P s 

(&tm) 

P 8 /P w 

T s 

(K) 

n s 

(cm) 

^c,w 

(MW/«?) 

q r,w 
(MW/m 2 ) 

30 

0.1141E-03 

0.1770 

17.77 

11938 

0.8026 

12.601 

0.5307 

40 

0.416 IE -03 

0.6200 

16.70 

12606 

0.8689 

21.193 

3.5369 

50 

0 . 1860E-02 

2.3840 

15.55 

12628 

0.8778 

31.959 

9 .6454 

60 

0. 5086E-02 

4. 1380 

16.67 

8881 

0.8786 

24. 967 

1.0094 

61 

0.5333E-02 

4.0860 

16.64 

8630 

0.8866 

23.314 

0.8945 

70 

0. 535GE-02 

2.4490 

14.87 

7290 

1.0150 

10.687 

0.3173 

80 

0.3397E-02 

1.0850 

13.18 

6507 

1.1430 

4.881 

0.1143 

100 

0. 1343E-02 

0.3143 

11.77 

5791 

1.2758 

1.769 

0.0327 
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Table 13. (Concluded.) 


(b) 


TIME 

(•) 

h 

w 

(kJ/kg) 

h . 

(kJ/kg) 

St 

Re 

30 

2. 045E+03 

82. 571E+03 

0.1056 

0. 1451E+04 

40 

2.045E+03 

79.543E+03 

0.5154E-01 

0.4664E+04 

50 

2.045E+03 

68. 739E+03 

0.2175E-01 

0. 1636E+05 

60 

2.045E+03 

43.342E+03 

0.1261E-01 

0 .3902E+05 

61 

2.045E+03 

40. so ie + o: 

0.1233E-01 

0.4089E+05 

70 

2.045E+03 

24.443E+03 

0.1255E-01 

0.3898E+05 

80 

2. 045E+03 

17.144E+03 

0. 1595E-01 

0.2371E+05 

100 

2.045E+03 

12.612E+03 

0.24255E-01 

0.9115E+04 


Table 14* Stagnation results (sphere cone, * 0*2 m, T w * 2,000 K) : 
atmosphere - 90% N 2 ♦ 10% CU 4 . 


(a) 


TIME 


P 

8 

P 8 ^ P « T S 

n . 

^c,w **r,w 

(•) 

(kg /. 3 ) 

(atm) 

(K) 

(c.) 

(MW/ ^ ) (MW/m 2 ) 


Atmosphere - 90% N 2 ♦ 10% CHi,, Trajectory I 


40 

0. 1265E-3 

0.1165 

18.56 

7259 

0.8029 

5.247 

0.0177 

50 

0.3473E-3 

0.3112 

17.95 

7460 

0.8440 

8.126 

0.0713 

60 

0. 1022E-2 

0.8455 

17.15 

7544 

0.8880 

11.637 

0.3569 

70 

0.3197E-2 

2.1108 

15.74 

7342 

0.9619 

13.498 

2.0398 

78 

0. 5157E-2 

2.443 

14.08 

6733 

1.0642 

9.687 

3.7697 

90 

0.55I2E-2 

1.512 

12.10 

5269 

1.2120 

3.543 

3.8894 

00 

0.3944E-2 

0.7999 

12.10 

4127 

1.2208 

1.683 

1.5915 




Atmosphere 

- 90% N 2 

+ 10% CHi*, Trajectory I 


140 

0.3812E-3 

0.2942 

17.78 

6953 

0.8586 

6.437 

0.0862 

150 

0.4687E-3 

0.3449 

17.49 

6895 

0.8731 

6.523 

0.1190 

170 

0..6712E-3 

0.4355 

16.75 

6729 

0.9111 

6.242 

0.2044 

175 

0. 7052E-3. 

0.4413 

16.58 

6665 

0.9201 

6.061 

0.2206 

180 

0.7213E-3 

0.4350 

16.39 

6596 

0.9300 

5.747 

0.2287 

200 

0.7855E-3 

0.4081 

15.63 

6330 

1.9723 

4.589 

0.2728 
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Table 14. (Concluded. ) 


(b) 


TIME 

(a) 

b w 

(kJ/kg) 

h . 

(kJ/kg) 

St 

Be 


Atmosphere 

- 902 Ns ♦ 

102 CU 4 , Trajectory I 

40 

2. 877E+03 

48.397E+03 

0.914CE-01 

0.1221E+04 

SO 

2.877E+03 

47.175E+03 

0.5365E-01 

0.3250E+04 

60 

2.877E+03 

43.595E+03 

0.2949E-01 

0.9229E+04 

70 

2.877E+03 

38.804E+03 

0.1554E-01 

0.2736E+05 

78 

2.877E+03 

25.002F+03 

0.1173E-01 

0.4218E+05 

90 

2.877E+03 

14.347E+03 

0.1007E-01 

0 .4460E+05 

100 

2.877E+03 

10.440E+03 

0.1182E-01 

0.340® +05 


Atmosphere - 902 N? + 10% CH^, Trajectory II 


140 

2.877E+03 

40.565E+03 

0.4904E-01 

0.3576E+04 

150 

2. 877E+03 

38.609E+03 

0.4352E-01 

0.435® +04 

170 

2.877E+03 

34.094E+03 

0.3545E-01 

0.6097E+04 

175 

2. 877E+03 

32.835E+03 

0.347CE-01 

0.637® +04 

180 

2.877E+03 

31 .651E+03 

0.3411E-01 

0.6494E+04 

200 

2.877E+03 

27.239E+03 

0.3175E-01 

0.6941E+04 


Table 15. Stagnation result# (sphere cone, • 0.2 u, • 
atmosphere - 90Z N 2 + 10Z CH^, Trajectory 111. 


(a) 


TIME 

(s) 

kg/®* 

P . 

(atm) 

P e /P„ 

T . 

(K) 

n e 

(cm) 

^c,»r 
(MW/m 2 ) 

90 

0.6142E-03 

0.1917 

13.24 

5159 

1.1301 

1.661 

100 

0.1076E-02 

0.3164 

12.81 

5210 

1.1614 

1.907 

no 

0.234SE-02 

0.6168 

12.31 

4974 

1.19!3 

2.915 

120 

0.3859E-02 

0.8309 

12.16 

4297 

1.2120 

1.864 

129 

0.5318E-02 

0.8997 

11.59 

3771 

1.2780 

1.330 

140 

0.6079E-02 

0.7416 

10.51 

3232 

1 .4083 

0.643 

150 

0.5643E-02 

0.5286 

9.74 

2790 

1.5168 

0.258 

160 

0.4666E-02 

0.3564 

9.40 

2397 

1.5725 

0.085 

170 

0.3693E-02 

0.2433 

9.39 

2081 

1.5811 

0.001 


2,000 K) 


(MW/m 2 ) 

0.2419 

0.5142 

1.3768 

1.8469 

1.0832 

0.2475 

0.0349 

0.0027 

0.0001 


Table IS. (Concluded.) 


(b) 


TIME 

(•) 

h 

V 

(kJ/kg) 

h . 

(kJ/kg) 

SC 

Be 

90 

2. 877E+03 

16. 30 7E +03 

0.3417E-01 

0.5213E+04 

100 

2.877E+03 

15.366E+03 

0.2477E-01 

0.8998E+04 

110 

2.8774+03 

13.694E+03 

0.1590E-01 

0.1949E+05 

120 

2.877E+03 

11.116E+03 

0.1190E-01 

0.331 1E+05 

129 

2.877E+03 

8.627E+03 

0.9913E-02 

0.4579E+05 

140 

2.877E+03 

6.099E+03 

0.8725E-02 

0.5154E+05 

150 

2.877B+03 

4.564E+03 

0.8114E-02 

0.4779E+05 

160 

2.877E+03 

3.611E+03 

0.7989E-02 

0.4040E+05 

170 

2. 877E+03 

3.022E+03 

0. 7906E-03 

0.3311E+05 


Table 16 


Stagnation results (sphere cone, • 0.2 m, ■ 2,000 K) : 
atmosphere - 90Z N2 + 102 CH4, Trajectory IV. 


(a) 


TIME 


P . 

P,/P- 

T . 

°s 


q r,w 

(a) 

(kg/«3) 

(atm) 


(R) 

(cm) 

(MW/m*) 

(MW/m 2 : 

60 

0.2927E-03 

0.1689 

16.81 

6233 

0.9092 

3.458 

0.0704 

70 

0 . 7469E-03 

0.4199 

16.06 

6469 

0.9482 

5.197 

0.2430 

80 

0.1689E-02 

0.8356 

14.92 

6468 

1.0130 

6.096 

0.7787 

90 

0.3726E-02 

1.4447 

13.37 

6172 

1.1155 

5.753 

2.4460 

07 

0.5162E-02 

1.5660 

12.33 

5601 

1.1946 

4.219 

3.8516 

110 

0.5786E-02 

1.0828 

11.83 

4002 

1 .2498 

1.647 

1.7776 

120 

0.4621E-02 

0.6529 

11.09 

3431 

1.3387 

0.780 

0.4077 

130 

0.3414E-02 

0.3972 

10.53 

3089 

1.4118 

0.389 

0.1017 

150 

0.1473E-02 

0.1381 

10.03 

2692 

1.4921 

0.107 

0.0093 
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Table 16. (Concluded.) 


(b) 


TIME 

(•) 

h 

V 

(kJ/kg) 

h . 

(kJ/kg) 

SC 

Re 

60 

2.877E+03 

30.243E+03 

0. 542GE-01 

0. 270® +04 

70 

2.877E+03 

29.538E+03 

0.3331E-01 

0.6663E+04 

60 

2.8774+03 

26.032E+03 

0. 2113E-01 

0. 1443E+05 

90 

2.877E+03 

20.401E+03 

0.1341E-01 

0.3034E+05 

97 

2.877E+03 

15. 934E+03 

0.1072E-01 

0.4140E+05 

110 

2.877E+03 

9.632E+03 

0.9168E-02 

0.4950E+05 

120 

2.877E+03 

7 . 124E+03 

0. 9871E-02 

0.3974E+05 

130 

2.877E+03 

5.776E+03 

0. 1068E-01 

0.2932E+05 

150 

2.877E+03 

4.574E+03 

0. 1308E-01 

0. 1277E+05 


Table 17. Stagnation reaulta (aphere cone, » 0.2 a, ■ 2,000 K): 
ataoaphere - 90% N 2 ♦ 10X CHt,, Trajectory V. 


(a) 


TIME 

(a) 

h * 

(kg/a 3 ) 

30 

0.1142-03 

40 

0.4165E-03 

30 

0.1862-02 

60 

0.5092E-02 

61 

0.5332-02 

70 

0.5357E-02 

80 

0.3401E-02 

100 

0.1345E-02 

120 

0.3482-03 


p, 0,/P. 


(atm) 

0.1765 

16.68 

0.6177 

15.73 

2.3770 

14.93 

4.1260 

15.96 

4.0730 

15.81 

2.4360 

13.89 

1.0810 

12.63 

0.3159 

12.46 

0.1075 

12.67 


T . 

(K) 

(ca) 

11326 

0.8531 

11869 

0.9243 

11542 

0.9876 

8162 

0.9483 

7973 

0.9569 

6650 

1.0772 

5650 

1.1720 

4468 

1.1865 

3949 

1.1793 


^c,a **r,w 


(MW/a 2 ) 

(MW/a 2 ) 

11.458 

0.3785 

19.458 

2.2628 

31.222 

6.2089 

23.683 

4.0362 

21.928 

4.3607 

9.004 

3.8994 

3.642 

2.3153 

1.228 

0.6681 

0.724 

0.3268 


Table 17. (Concluded.) 


(b) 


TIME 

(e) 

h 

V 

(kJ/kg) 

h 

s 

(kJ/kg) 

St 

Re 

30 

2.877E+03 

82.095E+03 

0.9761-01 

0.1349E+04 

40 

2.877E+03 

79.105E+03 

0.4812E-01 

0.4394E+04 

50 

2.877E+03 

68. 1986*03 

0. 21666-01 

0. 1616E+05 

60 

2.877E+03 

42.896E+03 

0.1235E-01 

0 .4308E+05 

61 

2.877E+03 

40. 33 IE +03 

0.1199E-01 

0.4505E+05 

70 

2.877E+03 

24.00VE+03 

0.1121E-01 

0.4358E+05 

80 

2.877k ♦03 

16.682E+03 

0. 130CE-01 

0.2756E+05 

100 

2.877E+03 

12.170E+03 

0.1917E-01 

0.1158E+05 

120 

2.877E+03 

10.607E+03 

0.2299E-01 

0. 7620E+04 


2,000 a): 


Table 18. 


Stagnation results (sphere cone, ■ 0.2 a, 
atmosphere - 90% H 2 + 10% Cl^, Trajectory VI 


T » 

w 


(a) 


TIME 


P . 

P« /p - 

T . 

°s • 

q c.a 

q r,w 

(s) 

kg/ffl? 

(ata) 


(K) 

(ca) 

(MW/a 2 ) (MW/n 2 ) 

60 

0.4396E-03 

0.5542 

16.32 

11512 

0.8800 

16.170 

1. 1565 

70 

0.1019E- 

1.1605 

15.91 

11160 

0.8993 

20.085 

1.5266 

73 

0.1227E-,i 

1 3348 

15.96 

10799 

0.8994 

20.276 

1.2076 

77 

0.1483E-02 

1.5028 

16.24 

10099 

0.8904 

19.715 

0.6897 

81 

0. 1672E-02 

1.5648 

16.80 

9241 

0.8722 

18.288 

0.3535 

90 

0.1715E-02 

1.3340 

17.34 

8059 

0.8666 

14.078 

0.2340 

100 

0.134 IE -02 

0.8811 

17.20 

7422 

0.8843 

9.543 

0.1415 

110 

0 . 1 109E-02 

0.6433 

16.87 

7036 

0.9055 

7.171 

0.1001 

120 

0. 1012E-02 

0.5296 

16.47 

6835 

0.9286 

5.798 

0.0866 

130 

0.9926E-03 

0.4630 

16.00 

6664 

0.9557 

4.929 

0.0813 




1 

1 
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Table 18. (Concluded. ) 


(b) 


TIME 

(«) 

h 

w 

(kJ/kg) 

h 

8 

(kJ/kg) 

St 

Re 

60 

2. 173E+03 

67.396E+03 

0.4813E-01 

0.4097E+04 

70 

2.173E+03 

60.923E+03 

0.3007E-01 

0.8400E+04 

73 

2.173E+03 

58.343E+03 

0.2707E-01 

0.986 IE +04 

77 

2.173E+03 

54.056E+03 

0.2434E-01 

0.1167E+05 

81 

2. 173E+03 

49.834E+03 

0.2272E-01 

0.1327E+05 

90 

2.173E+03 

41.342E+03 

0.2282E-01 

0.1396E+05 

100 

2. 173E+03 

34.823E+03 

0.2576E-01 

0.1097E+05 

110 

2.173E+03 

30.710E+03 

0.2845E-01 

0.9019E+04 

120 

2. 173E+03 

27.84GE+03 

0. 2962E-01 

0.8154E+04 

130 

2.173E+03 

25.520E+03 

0.2975E-01 

0.7894E+04 


Table 19. Downstream results (sphere cone, ■ 0.2 m, T w * 2,000 K) 
atmosphere - 98X N2 + 2Z CH^, Trajectory VI, Z * 196.3 km. 


N (MU/ nr ) 


0 

22. 

1.1 

21. 

1.2 

20. 

1.3 

20. 

1.4 

18. 

0 

21. 

1.1 

22.i 

1.2 

21. 

1.3 

20. 1 

1.4 

19. 


Enthal 


NO SLIP 

0.9568 2177 58250 10799 2000 15.95 0.0295 
0.9834 2177 57820 10724 2000 15.95 0.0282 
0.9890 2177 56240 10440 2000 16.09 0.0279 
1.0040 . 2177 53560 9904 2000 16.40 0.0268 
1.0268 2177 50080 9161 2000 16.92 0.0249 


BODY AND SHOCK SLIP 
1.0001 2384 58250 10799 2149 15.95 0.0295 

0.9743 2383 57940 10745 2149 15.91 0.0295 

0.9963 2383 56410 10471 2157 16.04 0.0283 

1.0181 2395 53780 9950 2159 16.33 0.0274 

1.0455 2400 50370 9221 2167 16.83 0.0255 





Table 19. (Concluded.) 


(b) 



q 

n c f w 

N 

8 

Enthalpy (kJ/kg) 

T s 

T v 


St 

S/R„ 

(MW/ nr 2 ) 

(cm) 

vail 

shock 

(K) 

(K) 







NO SLIP 





0 

22.59 

0.9386 

2381 

58330 

10811 

2151 

15.93 

0.0301 

0.1 

21.62 

0.9616 

2388 

57720 

10708 

2152 

15.98 

0.0288 

0.2 

20.77 

0.9793 

2382 

56200 

10433 

2149 

16.09 

0.0278 

0.3 

20.21 

1.0060 

2381 

53580 

9910 

2150 

16.40 

0.0269 

0.4 

18.71 

1.0430 

2392 

50250 

9 197 

2152 

16.89 

0.0251 




BODY AND SHOCK SLIP 




0 

22.32 

0.9674 

2177 

58330 

10812 

2000 

15.92 

0.0298 

0.1 

21.84 

0.9572 

2177 

57840 

10782 

2000 

15.94 

0.0291 

0.2 

20.93 

0.9853 

2177 

56360 

10462 

2000 

16.04 

0.0279 

0.3 

20.04 

1.0162 

2177 

53760 

9946 

2000 

16.34 

0.0267 

0.4 

18.86 

1.055b 

2177 

50440 

9236 

2000 

16.82 

0.0251 
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Table 20. Downstream results (sphere cone, ■ 0.2 m, • 2,000 K) : 
atmosphere - 982 N 2 + 22 CH^ , Trajectory VI, 2 * 241.8 km, 
e - 0.0515. 


(a) 



^C,W 

N 

3 

Enthalpy (IcJ/kg) 

T s 

T 

w 

P 8 /p - 

St 

s / r h 

(MW/n?) 

(cm) 

wall 

shock 

(K) 

(K) 







NO SLIP 





0 

17.27 

0.9126 

2179 

67540 

11509 

2000 

16.33 

0.0514 

0.1 

16.59 

0.9357 

2179 

66860 

11441 

2000 

16.33 

0.0493 

0.2 

16.16 

0.9567 

2179 

65070 

11254 

2000 

16.33 

0.0481 

0.3 

15.44 

0.9909 

2179 

62240 

10931 

2000 

16.33 

0.0459 

0.4 

14.49 

1.0374 

2179 

58320 

10409 

2000 

16.46 

0.0431 




BODY AND SHOCK SLIP 




0 

18.35 

0.9470 

2564 

67550 

11509 

2290 

16.33 

0.0546 

0.1 

18.16 

0.9317 

2570 

67020 

11456 

2279 

16.29 

0.0540 

0.2 

17.70 

0.9643 

2567 

65290 

11276 

2389 

16.28 

0.0526 

0.3 

16.95 

1.0023 

2577 

62420 

10952 

2310 

16.29 

0.0504 

0.4 

16.26 

1.0515 

2590 

58540 

10440 

2328 

16.42 

0.0484 
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Table 20. (Concluded. ) 


(b) 



^C,W 

N 

8 

Enthalpy (kJ/kg) 

T s 

T 

w 

P 8 /p - 

St 

s/r n 

(MW/rf) 

(cm) 

wall 

shock 

(K) 

(K) 







NO SLIP 





0 

18.52 

0.9141 

2570 

67520 

11506 

2289 

16.33 

0.0511 

0.1 

17.82 

0.9371 

2570 

66880 

11443 

2277 

16.33 

0.0530 

0.2 

16.03 

0.9584 

2572 

65100 

11258 

2285 

16.33 

0.0480 

0.3 

15.44 

0.9931 

2577 

62280 

10936 

2286 

16.33 

0.0462 

0.4 

14.39 

1.0399 

2586 

58370 

10416 

2291 

16.46 

0.0431 




BODY AND SHOCK SLIP 




0 

17.03 

0.9467 

2179 

67510 

11506 

2000 

16.33 

0.0506 

0.1 

17.54 

0.9313 

2179 

67020 

11456 

2000 

16.29 

0.0522 

0.2 

16.73 

0.9638 

2179 

65290 

11276 

2000 

16.28 

0.0498 

0.3 

15.97 

1.0016 

2179 

62410 

10951 

2000 

16.29 

0.0475 

0.4 

15.01 

1.0507 

2179 

58530 

10438 

2000 

16.42 

0.0446 
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Table 21. Downstream results (sphere cone, • 0.2 m, ■ 2,000 K) : 
atmosphere - 98% N 2 + 2% CH^, Trajectory VI, Z ■ 402.6 km, 
e ■ 0.286. 


(a) 



^c,w 

N 

8 

Enthalpy (kJ/kg) 

T s 

T w 

P s /P« 

St 

s/k b 

(MW/n£ ) 

(cm) 

wall 

shock 

(K) 

(K) 







NO SLIP 





0 

4.47 

0.6685 

2195 

71540 

10080 

2000 

18.79 

0.360 

0.1 

4.38 

0.6783 

2195 

70750 

10040 

2000 

18.80 

0.353 

0.2 

4.27 

0.6921 

2195 

68880 

9931 

2000 

18.71 

0.344 

0.3 

4.07 

0.7145 

2195 

65650 

9723 

2000 

18.60 

0.328 

0.4 

3.8C 

0.7470 

2195 

61460 

9412 

2000 

18.45 

0.306 




BODY AND SHOCK SLIP 




0 

3.80 

0.7457 

5636 

69810 

9999 

3749 

19.18 

0.322 

0.1 

3.92 

0.7235 

5702 

69670 

9988 

3785 

19.08 

0.333 

0.2 

3.57 

0.7736 

5743 

67910 

9879 

3803 

19.00 

0.303 

0.3 

3.26 

0.8195 

5916 

65010 

9686 

3884 

18.87 

0.278 

0.4 

2.28 

0.8725 

6065 

60860 

9371 

3945 

18.79 

0.245 


138 


Table 21. (Concluded. ) 


(b) 



^c,w 

N 

8 

Enthalpy (kJ/kg) 

T . 

T 

w 

P. /p - 

sc 

S/R„ 

(MW/n?) 

(cm) 

wall. 

shock 

(K) 

(K) 







NO SLIP 





0 

3.98 

0.6702 

6147 

71490 

10086 

2994 

18.80 

0.340 

0.1 

3.95 

0.6764 

6144 

70840 

10048 

3989 

18.78 

0.337 

0.2 

3.83 

0.6897 

6164 

68680 

9919 

3997 

18.70 

0.327 

0.3 

3.62 

0.7123 

6215 

65210 

9692 

4012 

18.57 

0.309 

0.4 

3.43 

0.7554 

6177 

61470 

9413 

3962 

18.45 

0.292 




BODY AND SHOCK SLIP 




• 

0 

3.98 

0.7417 

2195 

69730 

9995 

2000 

19.21 

0.321 

0.1 

4.18 

0.7203 

2195 

69630 

9986 

2000 

19.09 

0.337 

0.2 

3.84 

0.7669 

2195 

67870 

9876 

2000 

19.02 

0.310 

0.3 

3.55 

0.8098 

2195 

64930 

9680 

2000 

18.89 

0.287 

0.4 

3.19 

0.8598 

2195 

60710 

9359 

2000 

18.81 

0.258 
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Table 22. Downstream results (sphere cone, - 0.2 m, T w ■ 2,000 K) : 



atmosphere - 

982 N 2 ♦ 

22 CHt,, 

Trajectory VI, Z 

- 465.1 

kffl; 


e ■ 

0.524. 



(a) 





**c,w 

N 

8 

Enthalpy (kJ/kg) 

T a 

T w 


SC 

s/r n 

(MW/m 2 ) 

(cm) 

wall 

shock 

(K) 

(K) 







NO SLIP 




0 

3.45 

0.6042 

2208 

71530 

9545 

2000 

19.74 

0.934 

0.1 

3.39 

0.6113 

2208 

70830 

9510 

2000 

19.70 

0.918 

0.2 

3.32 

0.6217 

2208 

68860 

9407 

2000 

19.62 

0.896 

0.3 

3.17 

0.6396 

2208 

65650 

9226 

2000 

19.46 

0.857 

0.4 

2.96 

0.6654 

2208 

61320 

8947 

2000 

19.27 

0.801 




BODY AND SHOCK SLIP 




0 

1.76 

0.6294 

8970 

53530 

8347 

4532 

23.72 

0.527 

0.1 

2.05 

0.6039 

9955 

58030 

8768 

4647 

23.53 

0.624 

0.2 

1.72 

0.7023 

9583 

57220 

8697 

4600 

23.36 

0.519 

0.3 

1.55 

0.7676 

9225 

55800 

8556 

4551 

22.98 

0.465 

0.4 

1.35 

0.8398 

8631 

53000 

8249 

4462 

22.76 

0.401 


Table 22. (Concluded. ) 


(b) 



q c.w 

N 

8 

Enthalpy (kJ/kg) 

T 

s 

T 

w 


St 

s/ r n 

(MW/a? ) 

(cm) 

wall 

shock 

00 

0C) 







NO 8LIP 





0 

2.64 

0.6162 

12760 

71530 

9545 

4877 

19.74 

0.841 

0.1 

2.60 

0.6239 

12720 

70830 

9510 

4868 

19.70 

0.827 

0.2 

2.52 

0.6353 

12740 

68870 

9408 

4864 

19.62 

0.803 

0.3 

2.38 

0.6545 

12820 

65680 

9229 

4683 

19.46 

0.761 




BODY 

AND SHOCK 

SLIP 




0 

2.30 

0.6458 

2208 

55050 

8508 

2000 

24.97 

0.621 

0.1 

2.61 

0.5882 

2208 

56450 

8636 

2000 

24.18 

0.705 

0.2 

2.33 

0.6824 

2208 

56000 

8588 

2000 

23.87 

0.628 

0.3 

2.07 

0.7411 

2208 

54470 

8428 

2000 

23.58 

0.561 

0.4 

1.80 

0.8150 

2208 

52520 

8191 

2000 

23.06 

0.485 


Table 23. Downstream results (sphere cone, Rg - 0.2 m, T w ■ 2,000 K): 
atmosphere - 98% N 2 ♦ 2% CH 4 , Trajectory VI, Z - 497.6 km, 
e - 0.719. 




N 

s 

Enthalpy (kJ/kg) 

T 

s 

T 

w 

P # /P. 

St 

s/r n 

(MW/i ^ ) 

(cm) 

wall 

shock 

(K) 

(K) 







NO SLIP 





0 

22.49 

0.9386 

2381 

58330 

10811 

2151 

15.93 

0.0301 

0.1 

21.62 

0.9616 

2388 

57720 

10708 

2152 

15.98 

0.0288 

0.2 

20.77 

0.9793 

2383 

56200 

10433 

2149 

16.09 

0.0278 

0.3 

20.21 

1.0060 

2381 

53580 

9910 

2150 

16.40 

0.0269 

0.4 

18.71 

1.0430 

2392 

50250 

9197 

2152 

16.89 

0.0251 




BODY AMD SHOCK SLIP 




0 

22.32 

0.9674 

2177 

58330 

10812 

2000 

15.92 

0.0298 

0.1 

21.84 

0.9572 

2177 

57840 

10728 

2000 

15.94 

0.0291 

0.2 

20.93 

0.9853 

2177 

56360 

10462 

2000 

16.04 

0.0279 

0.3 

20.04 

1.0162 

2177 

53760 

9946 

2000 

16.34 

0.0267 

0.4 

18.86 

1.0556 

2177 

50440 

9236 

2000 

16.82 

0.0251 
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Table 24. Downstream results (sphere cone, ■ 0.2 m, • 2,000 K): 
atmosphere - 98X N 2 ♦ 2Z Cl^, Trajectory VI, Z ■ 530.8 km, 
e ■ 0.976. 



V 

N 

s 

Enthalpy (kJ/kg) 

T 

s 

T 

w 

P # /P. 

SC 

s/ s 

(MW/rf) 

(cm) 

wall 

shock 

00 

(K) 







NO 8LIP 





0 

3.10 

0.5704 

2240 

71490 

9032 

2000 

20.70 

2.908 

0.1 

3.05 

0.5755 

2240 

70760 

9000 

2000 

20.67 

2.859 

0.2 

2.97 

0.5849 

2240 

68770 

8909 

2000 

20.56 

2.792 

0.3 

2.85 

0.6004 

2240 

65590 

8752 

2000 

20.37 

2 675 

0.4 

2.67 

0.6237 

2240 

61250 

8508 

2000 

20.11 

2.507 




BODY 

AND SHOCK 

SLIP 




0 

2.19 

0.5937 

31590 

71490 

9032 

5383 

20.70 

3.561 

0.1 

2.13 

0.6007 

31470 

70770 

9001 

5373 

20.67 

3.448 

0.2 

1.98 

0.6167 

31360 

68880 

8914 

5366 

20.56 

3.212 

0.3 

1.77 

0.6435 

31220 

66060 

8776 

5344 

20.40 

2.855 

0.4 

1.49 

0.6767 

31000 

62300 

8571 

5322 

20.21 

2.392 
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